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Cf) ■ Abstract 

> : 

The hadronic charmless B decays into a scalar meson and a pseudoscalar meson are studied 
within the framework of QCD factorization. Based on the QCD sum rule method, we have derived 
\ the leading-twist light-cone distribution amplitudes of scalar mesons and their decay constants. Al- 

though the light scalar mesons /o(980) and ao(980) are widely perceived as primarily the four-quark 
bound states, in practice it is difficult to make quantitative predictions based on the four-quark pic- 
ture for light scalars. Hence, predictions are made in the 2-quark model for the scalar mesons. The 
short-distance approach suffices to explain the observed large rates of /o(980)K~ and /o(980)i^° 
that receive major penguin contributions from the b — > sss process. When /o(980) is assigned as 
a four-quark bound state, there exist extra diagrams contributing to B — > fo(980)K. Therefore, 
a priori the fo(980)K rate is not necessarily suppressed for a four-quark state /o(980). The pre- 
dicted B — > 0^(980)^ and clq (980)K~ rates exceed the current experimental limits, favoring a 
four-quark nature for ao(980). The penguin-dominated modes ao(980)K and ao(1450)K receive 
predominant weak annihilation contributions. There exists a two- fold experimental ambiguity in 
extracting the branching ratio of B~ — » K (1430)-7r~, which can be resolved by measuring other 
-Kq(1430)7t modes in conjunction with the isospin symmetry consideration. Large weak annihila- 
tion contributions are needed to explain the Kq(1A30)tt data. The decay B° — ► k + K~ provides a 
nice ground for testing the 4-quark and 2-quark nature of the k meson. It can proceed through 
VL-exchange and hence is quite suppressed if k is made of two quarks, while it receives a tree 
contribution if n is predominately a four-quark state. Hence, an observation of this channel at the 
level of > 10 may imply a four-quark assignment for the k. Mixing-induced CP asymmetries in 
penguin-dominated modes are studied and their deviations from sin 2/3 are found to be tiny. 
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I. INTRODUCTION 



The first charmless B decay into a scalar meson that has been observed is B — > fo(980)K. It 
was first measured by Belle in the charged B decays to K^tt^tt^ and a large branching fraction 
product for the fo(980)K ± final states was found 0] (updated in and 0) and subsequently 
confirmed by BaBar 0. Recently, BaBar has searched for the decays B — > cl^tt and B — » a$K for 
both charged and neutral ao mesons 0. Many measurements of B decays to other p-wave mesons 
such as ^(1430), / (1370), / ( 1500^ oi (1260), / 2 (1270), 00(1320) and K 2 *(1430) have also been 
reported recently by both BaBar 0,0, 0,0 and Belle P Tiol Illj|. The experimental results for 
the product of the branching ratios B(B — > SP) and B(S — ► P1P2) are summarized in Table [I] 
where S and P stand for scalar and pseudoscalar mesons, respectively. 

These measurements should provide information on the nature of the even-parity mesons. It 
is known that the identification of scalar mesons is difficult experimentally and the underlying 



structure of scalar mesons is not well established theoretically (for a review, see e.g. 
Studies of the mass spectrum of scalar mesons and their strong as well as electromagnetic decays 
suggest that the light scalars below or near 1 GeV form an SU(3) flavor nonet and are predominately 
the <?? states as originally advocated by Jaffe Q, while the scalar mesons above 1 GeV can 



be described as a qq nonet with a possible mixing with + qq and glueball states. It is hoped 
that through the study of B — > SP, old puzzles related to the internal structure and related 
parameters, e.g. the masses and widths, of light scalar mesons can receive new understanding. For 
example, it has been argued that a best candidate to distinguish the nature of the ao(980) scalar is 
B(B~ — > <2q 7t°) since the prediction for a four-quark model is one order of magnitude smaller than 
for the two quark assignment |l7j ]. 

One of the salient features of the scalar meson is that its decay constant is either zero or small of 
order — m u , m s — m^ iU . Therefore, when one of the pseudoscalar mesons in B — > PP decays is 
replaced by the corresponding scalar, the resulting decay pattern could be very different. Consider 
the decays B — > ao(980)7r as an example. It is expected that T(B + — ► ao~7r°) <C F(B + — ► a[)7r + ) 
and T(P° — ► Oq-k^) -C r(P° — ► Oq tt + ) as the factorizable contribution proportional to the decay 
constant of the scalar meson is suppressed relative to the one proportional to the pseudoscalar 
meson decay constant. This feature can be checked experimentally. 

Experimentally, BaBar 0] and Belle have adopted different approaches for parametrizing 
the non-resonant amplitudes in the 3-body decays B + — > K + ir + ir~ . Belle found two solutions 
with significantly different fractions of the B + — > Kq (1430)°7r + channel from the fit to K + it + it~ 
events. At first sight, it appears that the solution with the larger branching ratio, namely, B{B + — > 
Kq(1A30)°tt + ) ~ 45 x 10~ 6 [see Eq. (j2.4j) below], is preferable as it is consistent with the BaBar 
m _ent and supported by a pi—ological estimate h Q. However , since the counterpart 
of this dec ay i n the 2 pseudoscalar production, namely, B + — > K°t[ + has a branching ratio of order 
24 x 10" 6 hjf, one may wonder why the i^Q°7r + production is much more favorable than i^°7r + , 
while the i^Q 7r° mode is comparable to K°tt° (see Table HI)) ■ In this work, we shall examine the 
KqTt modes carefully within the framework of QCD factorization |2o| . 

Direct CP asymmetries in f$K and i^o(1430)7r modes have been measured recently by BaBar 



2 



and Belle (see Table III). Since direct CP violation is sensitive to the strong phases involved in the 
decay processes, the comparison between theory and experiment will provide information on the 
strong phases necessary for producing the measured direct CP asymmetries. 

The layout of the present paper is as follows. In Sec. II, we extract the absolute branching ratios 
of B — ► SP from the measured product of the branching ratios B(B — > SP) and B(S — ► PiPt)- 
The physical properties of the scalar mesons such as the quark contents, decay constants, form 
factors and their light-cone distribution amplitudes are discussed in Sec. III. We then apply QCD 
factorization in Sec. IV to calculate the branching ratios and CP asymmetries for B — * SP decays. 
Sec. V contains our conclusions. The factorizable amplitudes of various B — > SP decays are 
summarized in Appendix A. Based on the QCD sum rule method, the decay constants and the 
leading twist light-cone distribution amplitudes of the scalar mesons are evaluated in Appendices 
B and C respectively. 

II. EXPERIMENTAL STATUS 

The experimental results for the product of the branching ratios B{B — ► SP) and B{S — ► P1P2) 
are summarized in Tabled Here we shall try to determine B(B — > SP) given the information on 
B(S — > P\P2)- The absolute branching ratios for B — > fo(980)K and /n(980)7r depend critically on 
the branching fraction of /o(980) — > irir. For this purpose, we shall use the results from the most 
recent analysis of 21], namely, = 64 ± 8 MeV, T kR = 12 ± 1 MeV and r tot = 80 ± 10 MeV for 
/o(980). Therefore, 

B(/o(980) -» tt+tT) = 0.53 ±0.09, B(f (980) -> K + K~) = 0.08 ± 0.01 . (2.1) 

The obtained ratio r = 0(/ o (98O) -» 7r+7r-)/fi(/ (980) -» K + K~) ~ 7.1 is consistent with the 
result of r > 3.0l~Q'y inferred from the Belle measurements of B{B + — > /o(980)-ftT + — * it + it~K + ) 
and B(B+ -> /o(980)i^+ -> K+R-R+) (see Tabled. 

For an, we apply the Particle Data Group (PDG) average r(an — ► KK)/T(ao — > vrr/) = 0.183 ± 
0.024 Q to obtain 

S(a (980) -» r/vr) = 0.845 ±0.017. (2.2) 

Needless to say, it is of great importance to have more precise measurements of the branching 
fractions of fo and an- For ^(l^O) we have [22| 

B(K£° -»■ ir+Tr - ) = ^(0.93 ±0.10), fi(A" *° ^ K+^°) = ^(0.93 ± 0.10). (2.3) 

As noted in Tabled Belle found two solutions for the branching ratios of B + — ► i^g (1430)°7r + 
from the fit to 5+ - i^vr+vr- events Q. BaBar Q adopted a different approach to analyze 
the K + n + 7t~ data by parametrizing i^Q (1430)°7r + and the non-resonant component by a single 
amplitude suggested by the LASS collaboration to describe the scalar amplitude in elastic Kir 
scattering. As commented in jj], while this approach is experimentally motivated, the use of the 
LASS parametrization is limited to the elastic region of M(Kir) < 2.0 GeV, and an additional 
amplitude is still required for a satisfactory description of the data. Therefore, additional external 
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TABLE I: Experimental branching ratio products (in units of 10 6 ) of B decays to final states 
containing scalar mesons. The third error whenever occurred represents the model dependence. 



Mode 






BaBar [4, £, 1, 1, £, .9J Belle [2 , £ ljh JJ, .12] 


A Q"P 
JTV V C-L CI 


B(B + - 

LJ \ JJ 


— ± r?TC~^~ \ 
~ \j i\ j 




< 4.1 




< 4 1 


LJ \J-> 






9.3 ± 1.0 ± 0.5^07 a 


8.8 ± 0.8 ± 0.7+?! a 


O 1+O.8 
y • x 1.1 


B(B + - 

LJ \ J-J 




* ii ii / 




< 2.9 


<r 9 Q 

\ _ . C7 


B(B° - 

LJ I JJ 


-► fnCQRn 1 ) K°)B( fr,(QRQ) — ► 


7T^~ 7T ) 


5.5 ±0.7 ±0.6 


7.60 ± 1.66 ± 0.59+^ 


tj • iJ _l_ \J • O 


B(B + - 

LJ \ J-J 


-> fn(QR0)Tr+\B( fn(QRO) — > 


TT~^~ TT \ 
i\ ii J 


< 3.0 




/Id 


B(B + - 

LJ \ J-J 


-> n9.(QRCl')K + )B(nn( ( )R()) - 

* (in Wuu III J LJ \iJj[j \ U (J\J J 


-> r77r°l 


< 2.5 




\. _■ . (J 


B(B+ - 

LJ \ J-J 


~ IX q I i?OW 1 -i i. / LJ \ 11 [J I t/OU J 


— > A/7T ) 


< 3.9 






B(B + - 

LJ \ ±J 


-> ag(980)^ + )S(ao(980) - 




< 5.8 




^58 


ft( R° 

LJ \ -D — 


+ a (980)K+)B(a (980)- 


-» r/vr") 


< 2.1 


< 1.6 




B(B° - 


^ a§(980)^°)^(a (980) - 


>7?7r ) 


< 7.8 




< 7 S 


B(B° - 


^ a J (980)71^)0(00 (980)^ - 




< 5.1 


< 2.8 


< 2.8 


B(B+ - 


-> / (1370)A: + )B(/o(1370) 


— > 7T + 7r~) 


< 10.7 




< 10.7 


B(B+ - 


-> / (1370)7T+)e(/o(1370) 


— » 7T + 7T~) 


< 3.0 




< 3.0 


B(B+ - 


-> /o(1500)tf + )B(/ (1500) 


— > 7T + 7r~) 


< 4.4 




< 4.4 


B(B + - 


-> A * (1430)7r+),6(A * (1430) -> A+Tf-) 


34.4 ± 1.7 ± l.Sti\ h 


27.9 ± 1.8 ±2.6±|! c 












5.1±1.4±0.5lJ| c 




B(B° - 


^ A:^ + (1430)7r-)i3(A^ + (1430) -» A°vr+) 




30.8 ±2.4 ±2.4t^ 


30.8l|| 


B(B° - 


^ A * + (1430)O£(A* + (1430) -» ^ + ^°) 


11.2 ±1.5 ±3.5 




11.2 ±3.8 


B(B° - 


+■ A^°(1430)7r )e(A *°(1430) A+tt") 


7.9 ±1.5 ±2.7 




7.9±3.1 d 



a The previously published results are (9.2 ± 1.21^) x 10" 6 by BaBar A] and (7.6 ± 1.2+{;£) x 10" 6 by 







Belle 

^The BaBar result is for B+ -> (i^7r)^°Tr+ followed by (Kn)f -> K+tt" |7|. The (Kir)* ° component 
consists of a nonresonant effective range term plus the Aq°(1430) resonance itself. Using the knowledge of 
the composition of the -?£Tq (1430) component, BaBar obtained the branching ratio of B + — > _RTQ°(1430)7r + 
as shown in Eq. I|2.4|l . 

c Two solutions with significantly different branching ratios of the B + — > Kq (1430)°7r + channel but similar 
likelihood values were obtained by Belle from the fit to K + ir + ir~ events |2J. A new Belle measurement of 
K+it+tt- yields 32.0 ± 1.0 ± 2.4t^;^ for the larger solution Q. 

d The results B(B° -» X*+(1430)7r-)S(/i: o * + (1430) i^+vr ) = (5.1 ± 1.5±g;?) x lO" 6 and B{B° 
iq°(1430)7r )S(iq°(1430) -> K+tt-) = (6.1±i;|±g:|) x lO" 6 are quoted in ; 8| as Belle measurements. But 
they will not be included for the average as we cannot find these results in any Belle publications. 



information is needed in order to resolve the ambiguity in regard to the branching fraction of 
B + -» A(J(1430) 7r+: 

(37.0 ± 1.8 ± 1.9±°J ± 4.1) x 10~ 6 ; BaBar 



B(B + -» A^(1430)°7r+) = { (45.0 ± 2.9 ± 4.2± 1 |;^ ± 4.8) x 10" 

3.2 ± 2.2 ± 0.8±g;J ± 0.9) x 10~ 6 : 



Belle (solution I) [ 
Belle (solution II) 



(2.4) 
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TABLE II: Experimental branching ratios (in units of 10 6 ) of B decays to final states containing 
scalar mesons. 



Mode 




Br 


Mode 




Br 


B(B + - 
B(B+ - 
B{B+ - 
B{B + - 
B{B+ - 
B{B + - 
B{B+ - 


_|_ \ 

-> /o(980h+) 
-> / (980)^+) 
-> a[](980)7r + ) 
-> ag(980)i^+) 
-» a^(980)K°) 
-» ^°(1430)vr + ) 


< 4.1 

< 5.7 

1 7 1 +3-3 
i '- 1 -3.5 

< 6.9 

< 3.0 

< 4.6 
38.2±tl 


B(B U - 
B(B° - 
B(B° - 
B(B° - 
B(B° - 
B(B° - 


-> / (980)^ u ) 
■» 0^(980)^) 
- a°(980)^°) 
-» ao (980)K+) 
-» ^ (1430)vr ) 
K * + (1430)7r-) 


11.1 ± 2.4 
< 3.3 a 

< 9.2 

< 1.9 
12.7 ±5.4 

4 ' - z -6.9 



"Experimentally, one cannot separate B° from B decays, though theoretical calculations indicate T(B l 
a+n-) < T(B° -> a tt+) (see Table©. 



TABLE III: Experimental results of direct CP asymmetries in B decays to final states containing 
scalar mesons. 

Mode BaBar \7, .8, .9, 23] Belle [2, .3, 24, 25] Average 

B + -» / (980)A"+ 0.09 ± 0.10 ± 0.03lg;}o -0.077 ± 0.065 ± 0.030lg;^ -0.020l£$| 

#° -> /o(980)K° 0.24 ±0.31 ±0.15 -0.23 ± 0.23 ± 0.13 -0.06 ± 0.21 

# + -» /o(980)vr+ -0.50 ±0.54 ±0.06 -0.50 ± 0.54 

5+ -» if *°(1430)^+ —0.06 ± 0.03±g;g| 0.06 ±0.05^ -0.05±g;g| 

£° -> if *+(1430)7r- -0.07 ±0.12 ±0.08 -0.07 ±0.14 

B° -> i-q°(1430)7r -0.34 ±0.15 ±0.11 -0.34 ±0.19 



where the fourth error is due to the uncertainty on the branching fraction of .Kg (1430) — > Ktt [see 
Eq. (|2.3|) ]. For the BaBar result, the uncertainty on the proportion of the (Kit)^ component due 
to the i^o°(1430) resonance is also included in the fourth error. 

As shown in Sec. IV. B. 3, the aforementioned ambiguity can be resolved by measuring other 
i^o(1430)7r modes. The AI = penguin dominance implies, for example, the isospin relation 
T(B + — ► Kq°tt + ) = T(B° — > Kq + tt~). The recent measurements of the three-body decays B° — » 
K+vr+vr by BaBar [J and B° -> ^tt+tT by Belle Q yield 

, n x , N f (36.1 ±4.8 ±11.3 ±3.9) x 10~ 6 ; BaBaxB, , 

B (B° -» K'r (1430 +7T- = I) , 12 , ; - ' „ Hf 2.5 

v ov ; ; \ (49.7 ±3.8 ±3.8^|) x 10~ 6 ; Belle v ; 

It is clear that the isospin relation is well respected by both BaBar and Belle measurements of 
Kq°7t + and Kq + tt~ and that the smaller of the two solutions found by Belle (solution II) is ruled 
out. 

Experimental measurements of direct CP asymmetries for various B — > SP decays are shown 
in Table II I II We see that BaBar and Belle results for direct CP violation are consistent with zero. 
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III. PHYSICAL PROPERTIES OF SCALAR MESONS 



It is known that the underl y in g structure of scalar mesons is not well established theoretically 
(for a review, see e.g. 

QE1) It has been suggested that the light scalars below or near 
1 GeV-the isoscalars /o(600) (or a), /o(980), the isodoublet iQj;(800) (or k) and the isovector 
ao(980)-form an SU(3) flavor nonet, while scalar mesons above 1 GeV, namely,^ (1370), ao(1450), 
Kq(1430) and /o(1500)//o(1710), form another nonet. A consistent picture |l5( provided by the 
data suggests that the scalar meson states above 1 GeV can be identified as a conventional qq 
nonet with some possible glue content, whereas the light scalar mesons below or near 1 GeV form 
predominately a qqqq nonet 0, [3] with a possible mixing with + qq and glueball states. This 
is understandable because in the qq quark model, the + meson has a unit of orbital angular 
momentum and hence it should have a higher mass above 1 GeV. On the contrary, four quarks 
q 2 q 2 can form a + meson without introducing a unit of orbital angular momentum. Moreover, 
color and spin dependent interactions favor a flavor nonet configuration with attraction between 
the qq and qq pairs. Therefore, the + q 2 q 2 nonet has a mass near or below 1 GeV. This four-quark 
scenario explains naturally the mass degeneracy of /o(980) and ao(980), the broader decay widths 
of cr(600) and k(800) than / (980) and a (980), and the large coupling of / (980) and a (980) to 
KK. The four-quark flavor wave functions of light scalar mesons are symbolically given by 



80) 





a = uudd, /o = ss(uu + dd) / \/2, 

a{] = —j=(uu — dd)ss, clq = udss, = duss, 

k + = usdd, k° = dsuu, k° = sduu, k~ = sudd. (3-1) 

n 

This is supported by a lattice calculation \'2(i\ . 

While the above-mentioned four-quark assignment of light scalar mesons is certainly plausible 
when the light scalar meson is produced in low-energy reactions, one may wonder if the energetic 
/o(980) produced in B decays is dominated by the four-quark configuration as it requires to pick 
up two energetic quark-antiquark pairs to form a fast-moving light four-quark scalar meson. The 
Fock states of /o(980) consist of qq, q 2 q 2 , qqg, • • •, etc. Naively, it is expected that the distribution 
amplitude of /o(980) would be smaller in the four-quark model than in the two-quark picture. 

In the naive 2-quark model, the flavor wave functions of the light scalars read 

a = -^=(uu + dd), f = ss , 

a[j = —=(uu — dd), a<j = ud, clq = du, (3-2) 
v2 

K + = US, K° = ds, K° = sd, K~ = SU, 

where the ideal mixing for /o and a is assumed as /o(980) is the heaviest and a is the lightest one 
in the light scalar nonet. In this picture, /o(980) is purely an ss state and this is supported by the 
data of Df — > foir + and (p — > /o7 implying the copious /o(980) production via its ss component. 
However, there also exist some experimental evidences indicating that £$(980) is not purely an 
ss state. First, the observation of T(J/ip — > Jquj) ^T(J/^j — > /q^) .22}] clearly indicates the 
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existence of the non-strange and strange quark content in /o(980). Second, the fact that /o(980) 
and ao(980) have similar widths and that the /o width is dominated by tttt also suggests the 
composition of uu and dd pairs in /q(980); that is, /o(980) — ► tttt should not be OZI suppressed 
relative to ao(980) — * irrj. Therefore, isoscalars <r(600) and /o must have a mixing 

|/o(980)) = \ss) cos6 + \nn) Bin 6, k(600)) = -\ss) sin 9 + \nn) cos 9, (3.3) 

with nn = (uu + dd) /v2- 

Experimental implications for the fo — cr mixing angle have been discussed in detail in 

J/ip -► f <f>, / w[26] => 9 = (34 ± 6)° or = (146 ± 6)°, 

R = 4.03 ±0.14 [26] => 9 = (25.1 ± 0.5)° or = (164.3 ± 0.2)°, 
#= 1.63 ± 0.46 [26] = (42.3±jj;f) or # = (158 ±2)°, 

$ - /o7, /o - 77 [27] = (5±5)° or = (138±6)°, 

QCD sum rules and / data [28] 9 = (27 ± 13)° or 6» = (153 ± 13)°, 

QCD sum rules and a data [28] => 6» = (41 ± 11)° or 6» = (139 ± 11)°, (3.4) 

where R = 9 2 f oK + K - / g"f on + n - measures the ratio of the /o(980) coupling to K + K~ and -k + it~ . In 
short, 9 lies in the ranges of 25° < 9 < 40° and 140° < 9 < 165°. Note that the phenomenological 
analysis of the radiative decays 4> — > fo (980)7 an d /o(980) — > 77 favors the second solution, namely, 
9 = (138±6)°. The fact that phenomenologically there does not exist a unique mixing angle solution 
may already indicate that /o(980) and a are not purely qq bound states. 

Likewise, in the four-quark scenario for light scalar mesons, one can also define a similar fo — a 
mixing angle 

I /o (980)) = I nnss) cos <p + \ uudd) sin 4>, |<r(600)) = — \nnss) sincp + \uudd) cos <fi. (3.5) 
It has been shown that <ft = 174.6° [3^ . 



A. Decay constants 



To proceed we first discuss the decay constants of the pseudoscalar meson P and the scalar 
meson S defined by 

(-P(p)l<?27^75<?i|0) = -ifpPn, (S(p)\q2Jij,qi\0) = fsP^, (S\q 2 qi\0) = m s f s . (3.6) 

If the scalar meson is a four-quark bound state, it is pertinent to consider the interpolating current 
js, for example, 

jfo = ^^abc^dec[{ulC^S b ){u d ^Csl) + (u -> (£)], (3.7) 

with a,b,c - ■ ■ being the color indices and C the charge conjugation matrix. The coupling of the 
scalar meson S to the scalar current js is parametrized in terms of the scalar decay constant F$ 
defined by 

(S\j s \0) = V2F s m s . (3.8) 
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The neutral scalar mesons a, fo and a[] cannot be produced via the vector current owing to 
charge conjugation invar iance or conservation of vector current: 

U = ffo = f a o = 0. (3.9) 
For other scalar mesons, the vector decay constant fs and the scale-dependent scalar decay constant 
fs are related by equations of motion 

Vsfs = fs, with fi s = r -, (3.10) 

m 2 {fi) - mi(/i) 

where and m\ are the running current quark masses. Therefore, contrary to the case of pseu- 
doscalar mesons, the vector decay constant of the scalar meson, namely, fs, vanishes in the SU(3) 
or isospin limit. For example, the vector decay constant of Kq + (aj") is proportional to the mass 
difference between the constituent s (d) and u quarks; that is, the decay constants of Kq(1430) 
and the charged ao(980) are suppressed. In short, the vector decay constants of scalar mesons are 
either zero or small. 

For light scalar mesons, only two estimates of Fs in the four-quark scenario are available in the 
literature |3jJ, |32j and all other decay constant calculations are done in the 2-quark picture for light 



scalars. The results of Fs are 32 ] 



F a = (7.5 ± 1.0) MeV, F K = (1.6 ± 0.3) MeV, F fo = F ao = (1.1 ± 0.1) MeV. (3.11) 

We now turn to the model calculations in which the light scalar is assumed to be a two-quark 
bound state. Based on the finite-energy sum rule, Maltman obtained [3^] 

/ ao(980) = 1.1 ±0.2 MeV, / ao( i45o) =0.7±0.1MeV, f K . = 42 ± 2 MeV, (3.12) 

in accordance with the ranges estimated by Narison [3| 

/a (980) = 0.7 ~ 2.5 MeV, f K * = 33 ~ 46 MeV. (3.13) 

A different calculation of the scalar meson decay constants based on the generalized NLJ model 
yields 

/a (980) = 1-6 MeV, / ao{1450) = 0.4MeV, f K * = 31MeV. (3.14) 

Note that in ,33:] and [3| the ao decay constant is defined with an extra factor of (m s — m u ) / (m^ — 
m u ). We have taken the quark masses m s = 119 MeV, = 6.3 MeV and m u = 3.5 MeV at \i = 1 
GeV to convert it into our convention. Based on the QCD sum rule method, a recent estimate of the 
Kq (1430) scalar decay constant yields fx* = 427 ± 85 MeV at \i ~ 1 GeV j3] which corresponds 
to f K * = 34 ± 7 MeV. 1 

Because of the fo — a mixing, we shall treat /o and a separately. Just like the case of 77 and 77', 
each meson is described by four decay constants: 2 

(/o|««|0) = -j=m /o f%, (fo\ss\0) = m h f%, 

{a\uu\0) = ^=™*f2, <o-|»«|0> = mj°, (3.15) 



1 The estimate by Chernyak |l8| . namely, fx* — (70 ± 10) MeV, seems to be too large. 

2 Note that (ao|ss|0) = even when ao is a four-quark bound state. This is because ss is an isospin singlet 
while ao is an isospin triplet. 
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or 



</ n IHo> 



V2 

1 



m foff i 



{a n \uu\0) = ^rn a fS, 



(mss\0)=m fo f s fQ , 
(a s \ss\0) = m a f s a , 



(3.16) 



where /o\ c n = nn and f§, a s = ss. It follows that 

f% = f% S m9, cos 0, # = #coe0, /« = -/«sine. (3.17) 

Using the QCD sum-rule method, the scalar decay constant /| Q defined in Eq. (|3.16j) has been 



estimated in |38|] and |39| with similar results, namely, fj ~ 180 MeV at a typical hadronic scale. 
Taking into account the scale dependence of ff and radiative corrections to the quark loops in the 
OPE series, we have made a careful evaluation of the scalar decay constant in |40j using the sum 
rule approach. Our updated results for fj Q and f aQ of order 370 MeV at fx = 1 GeV (see Appendix 
B) are much larger than previous estimates. 3 Note that taking f ao = 1.1 MeV from Eq. ()3.12j) 
leads to f ao (lGeV) ~ 385 MeV, which is also very similar to our estimate. Therefore, a typical 
scalar decay constant of the scalar meson is above 300 MeV. In Appendix B we give a complete 
summary on the sum rule estimates of scalar meson decay constants. 



B. Light-Cone Distribution Amplitudes 

The twist-2 light-cone distribution amplitude (LCDA) <&s( x ) an d twist-3 &g{x) and <&g(x) for 
the scalar meson S made of the quarks (/29i are given by 



<5(p)|« 2 (^)gi(zi)|0) = m s / dxjte-*>+*r^$' s (x), 



{S{p)\q2{z2)a^qi{zi 



o 

Jo 



6 



with z = Z2 — zi, x = 1 — x, and their normalizations are 

dx® s (x) = fs, I dx$ s s (x) = [ dx$ s (x) = f s . 
io Jo Jo 

The definitions of LCDAs given in Eq. (|3,18|) can be combined into a single matrix element 



(3.18) 
(3.19) 



{S{p)\q2p{z2)qi a {zi 



X - £ dxe^ xp - Z2+ ^L^ s (x) + m s U s s (x) - <r^p 



a/3 



In general, the twist-2 light-cone distribution amplitude $5 has the form 



$s(x,/x) = f s (jJk)6x(l-x) 



B^)+Y^BMCU 2 {2x-l) 



m=l 



(3.20) 



(3.21) 



3 The decay constants fi and have been determined separately in ^jj] using the sum rule approach and 
they are found to be very close. Hence, for simplicity, we shall assume /| o = fj in the present work. 
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3 /2 

where B m are Gegenbauer moments and Cm are the Gegenbauer polynomials. The normalization 
condition (|3.19|) indicates 

B = /i s \ (3.22) 

where we have applied Eq. (|3.1Uf) and neglected the contributions from the even Gegenbaur mo- 
ments. It is clear that the Bq term is either zero or small of order uid — m u or m s — md tU , so are 
other even Gegenbaur moments [see also Eq. (jC3|) ]. For the neutral scalar mesons /o, a[| and a, 
Bq = and only odd Gegenbauer polynomials contribute. The LCDA also can be recast to the 
form 



®s(x,(x) = /s6x(l - x) 



(3.23) 



m=l 

which we shall use for later purposes. Since fis 3> 1 and even Gegenbauer coefficients are suppressed, 
it is clear that the LCDA of the scalar meson is dominated by the odd Gegenabuer moments. In 
contrast, the odd Gegenbauer moments vanish for the tt and p mesons. 

When the three-particle contributions are neglected, the twist-3 two-particle distribution am- 
plitudes are determined by the equations of motion, leading to 

(l-2x)^' s (x) = ^^-, (3.24) 

where use of Eq. (|3.18|) has been made. This means that we shall take the asymptotic forms 

*'s(x) = fs, = fs 6x(l - x), (3.25) 

recalling that it has been shown to the leading conformal expansion, the asymptotic forms of 
the twist-3 distribution amplitudes are the same as that for the pseudoscalar mesons j^. The 
corresponding light-cone projection operator of Eq. (|3.2U|) in momentum space can be obtained by 
assigning momenta 



20] 



k$ = xp^ + hf\ + -^—p^, K = x P v -kK + — p^, (3.26) 
2xp-p 2xp-p 

to the quark and antiquark in the scalar meson, where p is a light-like vector whose 3-components 
point into the opposite direction of p. As stressed in |2o| . the collinear approximation for the parton 
momentum (e.g. k\ = xp and &2 = xp) can be taken only after the light-cone projection has been 
applied. The light-cone projection operator of the scalar meson in momentum space then reads 



M^ = Up*$s(x)+m S j^mx)) , (3-27) 

^ "' 'a/3 



where use of Eq. (|3.25|) has been made. By comparison, the longitudinal part of the projection 
operator for the vector meson is given by |42| 

(*/iV\ i fv ( ( \ m vfv hfa i \ 
(M„ U = ~— [0v(x) - — j^jr*^) 

where the definitions for the twist-3 function <& v (x) and the transverse decay constant fy can 
be found in 42]. Therefore, the hard-scattering kernels for SP mesons in the final state can be 
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obtained from those for VP by performing the replacements fy&y{x) — > i ^,q( x) and myfy& v (x) — ► 
—ims&s(x), recalling that the normalization for $y an d & v is given by 42] 

f 1 dx<f>y(x) = 1, f 1 dx$ v {x) = 0. (3.29) 
Jo Jo 

Just as the decay constants for /o(980) and a, their LCDAs should also be treated separately. 

The twist-2 and twist-3 distribution amplitudes and (q = n, s) 4 , respectively, are given 

by 

(S^(p)\n(z)^n(0)\0) = Pft [* dxe**-**P(x), 

Jo 

(S^(pMz)^s(0)\0) = P, f 1 dxe^ g ^\x), 



o 



(,S (n) (p)|n(z)n(0)|0) = C dxe ixp - z ^ )s (x), 

Jo 

{S^(p)\s(z)s(0)\0) = m { s s) [ dxe ixp - z <$>f s {x). (3.30) 



They satisfy the relations <&s( x ) = ~~ ^s(l ~ x ) due to charge conjugation invariance (that is, the 
distribution amplitude vanishes at x = 1/2) and $g(x) = $5(1 — x) so that 

( l dx<S>f s \x) = 0, ( l dx$f s) \x) = /"' s , (3.31) 
Jo Jo 

with fg' s being defined in Eq. ()3.16j) . Hence, the light-cone distribution amplitudes for S = fo,a 
read 

*P'\x,ri = j»>'6x(l-x) J2 Bt s) ^)CU\2x-l). (3.32) 

m=l,3,5,-" 

The LCDAs are 

* /o (a;,/i) = cos 6» + $^ } sin0, $ ff (x,n) = -<S>l s) sm9 + $^ n) cos 9. (3.33) 

Since the Bq term in the LCDA for the charged ao is of order — m u , it can be safely neglected. 
Hence, in practice we shall use the same LCDA for both neutral and charged ao scalar mesons. 

Based on the QCD sum rule technique, the Gegenbauer moments in Eq. (|3.32|) have been 
evaluated in [40J up to m = 5. For an updated analysis, see Appendix C. Note that our result 
f ao Bi° = —340 MeV is much larger than the estimate of \fBi\ m 100 MeV at \i = m& inferred 
from the analysis in |3] (see Eq. (52) of 

For pseudoscalar mesons, the asymptotic forms for twist-2 and twist-3 distribution amplitudes 
for pseuodscalar mesons are 

$ P (x) = f P 6x(l-x), $ p p (x) = fp, $$(x) = f P 6x(l - x). (3.34) 



4 The quark flavor s should not be confused with the superscript s for the twist-3 LCDA <I> S (x) . 
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TABLE IV: Form factors of B — > tt, K, ao(1450), Kq (1430) transitions obtained in the covariant 



light-front model 45] 



7 ' 

b 


F(0) 


pi/ 2 \ 
" WmaxJ 


a 


7 

b 


b 


F(0) 




a 


7 




pB-K 


0.25 


1.16 


1.73 


0.95 


Ft 


0.25 


0.86 


0.84 


0.10 


pBK 
j-,Ba (U50) 


0.35 
0.26 


2.17 
0.68 


1.58 
1.57 


0.68 
0.70 


F?« 

^0(1450) 


0.35 
0.26 


0.80 
0.35 


0.71 
0.55 


0.04 
0.03 


F BK* 


0.21 a 
0.26 


0.52" 
0.70 


1.66" 
1.52 


1.00" 
0.64 




0.21" 
0.26 


0.33" 
0.33 


0.73" 
0.44 


0.09" 
0.05 




0.21" 


0.52" 


1.59" 


0.91" 




0.21" 


0.30" 


0.59" 


0.09" 



"Form factors obtained by considering the scalar meson above 1 GeV as the first excited state of the 
corresponding light scalar meson. 



C. Form factors 



Form factors for B — > P, S transitions are defined by 



H 



(P(p')\V,\B(p)) 
(S(p')\A,\B(p)) 



B 



m. 



% ) F? p tf) + 



m 



B 



ra 



P„ 



m 



B 



rrir 



(3.35) 



where P„ 



(p — p')^. As shown in 0|, a factor of (— i) is needed in B — > S 
-> S form factors to be positive. This also can be checked from heavy 



(p + p')m> % = 

transition in order for the B 
quark symmetry j^]. 

Various form factors for B — > S transitions have been evaluated in the relativistic covariant 
light-front quark model j^fj. In this model form factors are first calculated in the spacelike region 
and their momentum dependence is fitted to a 3-parameter form 

F(0) 



F(q 2 



1 — a(q 2 /m 2 B ) + b(q 2 /m 2 B ) 2 



(3.36) 



The parameters a, b and -F(O) are first determined in the spacelike region. This parametrization is 
then analytically continued to the timelike region to determine the physical form factors at q 2 > 0. 
The results relevant for our purposes are summarized in Table IIVI Note that the calculation of 
B to scalar meson form factors in coauthored by two of us is for the case where the scalar 
meson is made of qq' quarks. Since it is possible that K q (1430), a® (1450), /o (1500) are the first 
excited states of k, ao(980) and /o(980), respectively, we also extend the calculation to the case 
where -fCg(1430) and ao(1450) are first excited states by working out their wave functions from a 
simple-harmonic-oscillator-type potential. The resultant form factors are shown in Table fTVl 

Assuming that the light scalar mesons are the bound states of qq, form factors for B to light 
scalar mesons also can be estimated in this approach. Taking the decay constants of /o(980) and 
ao(980) estimated in Appendix B, it is found that the form factor of B to /o(980) or ao(980) is 
of order 0.25 at q 2 = 0. Therefore, the form factor ^?^ a °( 980 ) j s no ^ necessarily smaller than F BlT . 
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This is understandable because the ao(980) distribution amplitude peaks at x ~ 0.25 and x ~ 0.75 
while the pion LCDA peaks at x = 1/2. As pointed out in [2], since <J> ao is more pronounced 
towards the endpoints x = and x = 1, it can have a greater overlap with the highly asymmetric 
wave function of the B meson than the pion wave function can. Consequently, the B to ao(980) 
transition form factor is anticipated to be at least of the same order as the B — > tt case. Note that 
based on the light-cone sum rules, Chernyak 18] has estimated the B — > ao(1450) transition form 
factor and obtained i?^ a °( 1450 ) (o) = 0.46, while our result is 0.26 and is similar to the B — > it form 
factor at q 2 = 0. We will make a comment on this when discussing the decay B° — > o,|(1450)7r~ in 
Sec. IV.B. 



IV. B -> SP DECAYS 



A. Decay amplitudes in QCD factorization 



We shall use the QCD factorization approach 20, 42] to study the short-distance contributions 
to the decays B -> f (980)K, ETq(1430)7t, and a 7r, a K for a = a (980) and a (1450). In 
QCD factorization, the factorization amplitudes of above-mentioned decays are summarized in 
Appendix A. The effective parameters a? with p = u, c in Eq. I)A6|) can be calculated in the QCD 
factorization approach [2(3]. They are basically the Wilson coefficients in conjunction with short- 
distance nonfactorizable corrections such as vertex corrections and hard spectator interactions. In 
general, they have the expressions 



3,0 



a p i {M 1 M 2 ) 



. Cj±l C;±i CFa s i 



N r 4vr 



Air 2 

Vi(M 2 ) + —H i {M l M 2 ) + Pf{M 2 



(4.1) 



where i = 1, • • • , 10, the upper (lower) signs apply when % is odd (even), c, are the Wilson coefficients, 
C F = (N 2 - 1)/(2N C ) with iV c = 3, M 2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi{M 2 ) account for vertex corrections, Hi(M\M 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and Pi{M 2 ) for penguin contractions. The vertex and penguin corrections 
for SP final states have the same expressions as those for PP states and can be found in 
Using the general LCDA 



20 



<&m{x,ij) = /m6x(1 



l + E«nV)^ 2 (2* 



l) 



71=1 



with a n = [J-sBn for the scalar meson [see Eq. (|3.23p ] and applying Eq. 
corrections, we obtain (apart from the decay constant /m) 



(4.2) 

(37) in for vertex 



Vi{M) 



12 In 



18 



for % = 1 - 4, 9, 10, 
Vi{M) = 



mh 

12 In — + 6 



3in + 



3iir 



11 

T 



n 

y 



3i7T 



ai 



21 
20 



af + 



3Z7T 



Oil 



21 
20 



a 2 



79 
36 



79 
36 



3 J 3 



^ af + ••• 



(4.3) 



3 J 



(4.4) 
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for i = 5, 7 and Vi{M2) = —6 for i = 6, 8 in the NDR scheme for 75. The expressions of Vi(M) up 



to the a 2 term are the same as that in 

As for the hard spectator function H, it reads 



H i (M l M 2 ) 



1 



/ M2 F SMl (0)m| 



JO P JO £ JO ?7 



(4.5) 



for t = 1 - 4, 9, 10, 
Hi(M x M 2 ) = - 



fM 2 F mh (0)m% 



1 dp 



$b{p) 



1 d£ 



<*>M 2 (6 



1 dr\ 
»j 



oP ' •/.'! i •/() '/I x £ 

(4.6) 

for £ = 5, 7 and Hi = for z = 6,8, where £ = 1 — £ and fj = 1 — 77, $m (^m) is the twist- 2 
(twist-3) light-cone distribution amplitude of the meson M. The ratios and are defined 



in Eqs. QA5|I and (|A2|) . As shown in Appendix A, the factorizable amplitudes ^4p5 and A$p have 
an opposite relative sign [see Eq. (|A4|) ] and one has to replace by — when Mi is a scalar 
meson. This amounts to changing the sign of the first term in the expression of Hi{M\M 2 ) for a 
scalar meson M\. 

Weak annihilation contributions are described by the terms bi, and &i,EW i n Eq. (|A6|) which 
have the expressions 

Cf 
N 2 

c 

c f \_ Ai , _ A f 



b 2 



Cf A i 



Cf 
Cf 

&4,EW = ^2 



b 3 = 
bi 



c 3 A\ + c 5 (4 + A{) + N c c 6 A f 3 



c±A\ + c e A 



>3,EW 



A\ + c 7 (4 + A{) + N c c 8 Al 



c 10 A\ + c 8 Al 



(4.7) 



where the subscripts 1,2,3 of Aj^f denote the annihilation amplitudes induced from (V — A)(V — A), 
(V — A)(V + A) and (S — P)(S + P) operators, respectively, and the superscripts i and / refer to 
gluon emission from the initial and final-state quarks, respectively. Their explicit expressions are 
given by 



a\ 

A' 

Af 
J ^3 

M 



$M 2 (x)$ Ml (y) 
-$M 2 (x)$M 1 (y) 



y(l-xy) 
1 

y(l-xy) 
1 



x'y 
1 

Wy 



x(l — xy) xy 2 

1 _|_ 1 

x(l— xy) xy 2 



+ rf 1 rf 2 * m2 (x)$ mi (y) 3 



f^M 2 (x)$ mi (y) 



2y 



x) - , 2x 

I xy(l-xy) 



2x 

xy{l-xy) ) ' 



for M X M 2 = PS, 
for MiM 2 = SP, 
for MiM 2 = PS, 
for MiM 2 = 5P, 

for M X M 2 = PS, 
for MiM 2 = SP, 



rfi$M 2 (x)$ mi (y) ^ - r^$ Ml (y)$ m2 (x) 2 -^f) ; for M,M 2 = PS, 



x ~m2\™J~mi\») x 2 y 'X 
.Mifls.. f„,>» 2(1+3;) 



rf 1 $A/ 2 (x)$ mi (y) =^ - r x 



¥ 2 ®MAy)® 



in 2 



[x)^-); for M 1 M 2 = SP, 



A' 



(4.8) 
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where f ■ ■ ■ = ira s Jq dxdy, x = 1 — x and y = 1 — y. Note that we have adopted the same convention 
as in |42| that Mi contains an antiquark from the weak vertex with longitudinal fraction y, while 
M2 contains a quark from the weak vertex with momentum fraction x. 

Using the asymptotic distribution amplitudes for pseudoscalar mesons and keeping the LCDA 
of the scalar meson to the third Gegenbaur polynomial in Eq. (|3.23|) . the annihilation contributions 
can be simplified to 



A\(PS) 
Ah(PS) 

Al(PS) 

A f 3 (PS) t 
for M X M 2 



+ 




Bt(3X A + 4 - vr 2 ) + £3 ( WX a + 



23 10 , 

7T' 

18 3 
2956 100 



r S r P y2 



B X (X A + 29 - 3tt z ) +B 3 [X a + — 5- 



7T 



AX A + 4 + — 



9 



+ IOS3 Xj 



r s r p Xl 



19 v 191 vr 2 



Qfpfs^ot s X A W x us 
PS, and 

A\{SP) = Al(PS) 
Al(SP) = -Ai(PS) 



Bi(6X A -ll) + S 3 (20X A 



187\ 



r s y {2X A - r 



(4.9) 



Ai(SP) = A\(PS), 
4(SP)=4(PS), 



(4.10) 



for M1M2 = SP, where the endpoint divergence X A is defined in Eq. (|4.11|) . As noticed in passing, 
for neutral scalars a, fa and a§, one needs to express fs r x by fs^x an< ^ fsl^s by fs- Numerically, 
the dominant annihilation contribution arises from the factorizable penguin-induced annihilation 
characterized by A^. Physically, this is because the penguin- induced annihilation contribution is 
not subject to helicity suppression. 

Although the parameters ai{i 7^ 6, 8) and a^^r x are formally renormalization scale and 75 scheme 
independent, in practice there exists some residual scale dependence in a.;(/i) to finite order. To be 
specific, we shall evaluate the vertex corrections to the decay amplitude at the scale fi = In 
contrast, as stressed in [2(3], the hard spectator and annihilation contributions should be evaluated 
at the hard-collinear scale fj,h = y/JjXh with A^ ~ 500 MeV. There is one more serious complication 
about these contributions; that is, while QCD factorization predictions are model independent 
in the mj — * 00 limit, power corrections always involve troublesome endpoint divergences. For 
example, the annihilation amplitude has endpoint divergences even at twist-2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and posses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow 20] to parameterize the endpoint divergence Xa = Jq dx/x in the annihilation 
diagram as 



X A = \n 



( mB \ 

UJ 



(4.11) 



with the unknown real parameters pa and (ft a- Likewise, the endpoint divergence Xjj in the hard 
spectator contributions can be parameterized in a similar manner. 
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TABLE V: Branching ratios (in units of 1CP 6 ) of B decays to final states containing scalar mesons. 
The theoretical errors correspond to the uncertainties due to (i) the Gegenbauer moments -61,3, the 
scalar meson decay constants, (ii) the heavy-to-light form factors and the strange quark mass, and 
(hi) the power corrections due to weak annihilation and hard spectator interactions, respectively. 
The predicted branching ratios of B — ► /o(980)-fr, /o(980)-7r are for the f^ — v mixing angle = 155°. 
For light scalar mesons /o(980), ao(980) and k we have assumed the 2-quark content for them. The 
scalar mesons ao(1450) and -K"q(1450) are treated as the first excited states of ao(980) and k, 
respectively, corresponding to scenario 1 explained in Appendices B and C. Experimental results 



are taken from Table ITU 



"TV /f 1 

Mode 




1 neory 


Expt 


TV JT 1 

Mode 




I neory 


Expt 


B- -> 


fo(980)K- 


1 K R+0.3+4.7+5.4 
— 0.3—3.3 — 2.4 


17.1+H 

—3.5 


5°- 


fo(980)K U 


-,0 0+0.2+4. 1+4.5 

— 0.2—2.9—2.1 


11.2 ±2.4 


B -> 


/o(980)tt- 


n Q+0.0+0.3+0.2 
u - y -0.0-0.2-0.0 


< 5.7 




/o(980)vr° 


n n o+0. 01+0. 03+0. 08 
u - uo -0.01-0.00-0.01 




B~ -» 


0^(980)^" 


9 9+0.7+0.7+7.6 
z - z -0. 5-0.5-1. 7 


< 3.0 


B°^ 


a£(980)K- 


a q+1.3+1.4+14.8 
^•°-l. 1-1.0- 3.4 


< 1.9 


B~ -> 


a (980)iT° 


A q+1.4+1.8+16.1 
4 - y -l.l-1.2- 4.0 


< 4.6 


5°-> 


a%(980)K° 


9 ^+0.7+0.9+7.9 
z -^-0. 6-0.6-2.0 


< 9.2 


B~ -» 


ag(980)vr- 


a+0. 2+1.0+0.4 
°-^-0. 2-0.8-0. 4 


< 6.9 


Tf^ 


a(]"(980)7r- 


7 r+O. 7+2.0+2.2 
'• D -0. 6-1.8-1. 6 


< 3.3 a 


B~ -> 


a (980)vr° 


n 9+0.1+0.0+0.2 
u - z -o. 1-0.0-0.1 




B°^ 
U° — > 


0^(980)^+ 
a° (980)ir° 


n k+0.2+0.1+0.7 
u - D -0.1-0. 1-0.3 

n 9+0.1+0.0+0.1 

u - z -0.1-0.0-0.0 




B -» 


a%(U50)K- 


k fi +2.2+3. 5+8.6 
°- D -l. 7-1. 9-5.2 




Tf^ 


a^(U50)K~ 


H 1+4.4+6.9+17.1 
ii - i -3. 4-3.8-10. 2 




B~ -> 


a (U50)K° 


-i/i 1+5.0+8.2+18.9 
i ^- 1 -3.9-4.6-14.0 




B°^ 


a$(1450)K° 


fi fi +2. 3+3.9+9.0 
u - u -1.9-2.2-6.6 




B~ -> 


ag(1450)vr- 


4 1+0.5+1.1+1.3 
^-O. 4-1.0-1.1 




Tf^ 


0^(1450)71-- 


10 q+2.3+2.4+10.0 
lZ.y_ 2 0-2.2- 7.5 




B- -» 


a (1450)vr° 


n fi +0.2+0.1+0.5 
u - D -0.2-0. 1-0.3 




B°^ 
B°^ 


a|T (1450)tt + 
ag(1450)vr° 


n 1+0.1+0.1+0.9 
u - i -o.i-o. 0-0.0 
n 0+0.2+0. 1+0.2 
u -°-o.i-o. 1-0.1 




B~ -> 


Zo°(1430)vr- 


1 n +0.8+2.0+19.5 
1 ' u -0.5-0.7- 0.9 


38.2±ti 


B°^ 


i^0"(1430)7T+ 


1 1+0.8+2.1+17.7 
i - i -0. 5-0.9- 1.0 


47 9+ 5 - 6 
4 ' - z -6.9 




£T *-(1430)^° 


n Q+0.3+0.8+8.9 
u -°-0.2-0.2-0.3 




B°^ 


Zo°(1430)7r° 


n r+0.4+1.0+8.8 
u - D -0. 3-0.5-0. 5 


12.7 ±5.4 



°The cited upper limit 3.3 x 1CT 6 is for B -> (980)71-+. 



Besides the penguin and annihilation contributions formally of order 1/mt, there may exist 
other power corrections which unfortunately cannot be studied in a systematical way as they are 
nonperturbative in nature. The so-called "charming penguin" contribution is one of the long- 
distance effects that have been widely discussed. The importance of this nonpertrubative effect 
has also been conjectured to be justified in the context of soft-collinear effective theory [4f|. More 
recently, it has been shown that such an effect can be incorporated in final-state interactions 42]. 
However, in order to see the relevance of the charming penguin effect to B decays into scalar 
resonances, we need to await more data with better accuracy. 



B. Results and discussions 



While it is widely believed that /o(980) and ao(980) are predominately four-quark states, in 
practice it is difficult to make quantitative predictions on hadronic B — > SP decays based on the 
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TABLE VI: Same as Table IVl except that the mesons ao(1450) and Kq (1450) are treated as the 
lowest lying scalar states, corresponding to scenario 2 explained in Appendices B and C. 



Mode 



Theory Expt 



Mode 



Theory Expt 



w 

B° 
S° 
B° 
E° 
if 
B° 



a^(U50)K~ 0.2 



B~ -> 

B- Oq (1450)i? 
B- -> ag(1450)vr- 
B -> Oq (1450)tt° 



■n 



o (1430)vr- 
(1430)vr° 



0.2+0.1+17.6 
0.0-0.1- 0.0 
n 1+0.6+0.3+35.9 

u - i -o.o-o.o- 0.1 

9 c+0.3+0.9+1.2 
z -°-0. 3-0.7-0. 8 
1+0.4+0.1+1.1 



1.1" 



0.3-0.1-0.6 



n o+ 10 - 3 

iJ -- u - 6.0 



5.3 



7.5+49.9 no o- 
3.5-10.1 <J°- Z - 
+4.7+1.6+22.3 
2.8-1.7- 4.7 



•4.6 
•4.5 



• af(1450)A- 

• aj](1450)Z° 

• af(1450)7r- 
a(7(1450)vr+ 

• a§(1450)7T° 

Zo°(1430)7T 6.4 



1-0.5+0.3+36.4 
-0.0-0.1- 0.0 
1+0.3+0.2+17.7 
■ L -0.0-0.0- 0.0 
1-1.0+1.3+6.8 
-0.9-1.0-1.8 
^0.2+0.2+2. 6 
-0.2-0.2-0.3 
t-0.3+0.1+0. 5 
1 -0.2-0.0-0.3 
-q 3+9.4+3.7+45.8 



0.31 
0. 

3.r 

o.5i 
o.7: 



-5.8-3.7 
1-5. 4+2.2- 
-3.3-2.1 



9.9 
-26.1 
- 5.7 



47.2 



+5.6 
6.9 



12.7 ±5.4 



four-quark picture for light scalar mesons as it involves not only the unknown form factors and 
decay constants that are beyond the conventional quark model but also additional nonfactorizable 
contributions that are difficult to estimate (an example will be shown shortly below). Hence, we 
shall assume the two-quark scenario for /o(980) and ao(980). 

For form factors we shall use those derived in the covariant light-front quark model |45|. For 
CKM matrix elements we use the updated Wolfenstein parameters A = 0.825, A = 0.2262, p = 0.207 
and fj = 0.340 ^^|. For the running current quark masses we employ 

m h {m h ) = 4.2 GeV, m 6 (2.1 GeV) = 4.95 GeV, m 6 (l GeV) = 6.89 GeV, 

m c (m b ) = 1.3 GeV, m c (2.1 GeV) = 1.51 GeV, 

m s (2.1GeV) = 90MeV, m s (lGeV) = 119 MeV, 

m d (lGeV) = 6.3 MeV, m u (l GeV) = 3.5 MeV. (4.12) 
The strong coupling constants are given by 

a s (2.1 GeV) = 0.303, a s (l GeV) = 0.517, (4.13) 



221. 



corresponding to the world average a s (mz) = 0.1213 

The calculated results for branching ratios and CP asymmetries are exhibited in Tables fVl Villi 5 
In these tables we have included theoretical errors arising from the uncertainties in the Gegenbauer 
moments B\ ^ (cf. Appendix C), the scalar meson decay constant fs or fs (see Appendix B), the 
form factors F BP,BS , the quark masses and the power corrections from weak annihilation and hard 
spectator interactions characterized by the parameters Xa and Xh, respectively. For form factors 
we assign their uncertainties to be 5F BP ' BS (0) = ±0.03, for example, ^^(0) = 0.35 ± 0.03 and 
Fq K °(0) = 0.26 ± 0.03. The strange quark mass is taken to be m s (2GeV) = 90 ± 20 MeV. For 
the quantities Xa and Xh we adopt the form (|4.11|) with pa,H < 0.5 and arbitrary strong phases 
4>A,H- Note that the central values (or "default" results) correspond to pa,h = and 4>a,h = 0. 



5 B decays into light scalar mesons are not listed in Tables IVll and IVTTT1 as we do not have a handle for light 
scalars made of four quarks as explained in the text. 
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To obtain the errors shown in Tables IVIVIII1 we first scan randomly the points in the allowed 
ranges of the above six parameters in three separated groups: the first two, the second two and 
the last two, and then add errors in each group in quadrature. Therefore, the first theoretical 
error shown in the Tables is due to the variation of ^ and fs, the second error comes from the 
uncertainties of the form factors and the strange quark mass, while the third error from the power 
corrections due to weak annihilation and hard spectator interactions. 

Just like the B decays into PP or VP final states in the QCD factorization approach 2(], 4^ |. 



the theoretical errors are dominated by the l/mj, power corrections due to weak annihilation. 
However, it is clear from Tables IVIVII that the theoretical uncertainties in decay rates due to weak 
annihilation in some B — > SP decays, e.g. B — > 00(980)^, ao(1450)if and Kq(1430)it can be much 
larger than the "default" central values, while in B — > PP or VP decays, the errors due to Xah are 



comparable to or smaller than the central values (see e.g. Table 2 of 42j). This can be understood as 
follows. Consider the penguin-induced annihilation diagram for B — > PP. Its amplitude is helicity 
suppressed as the helicity of one of the final-state mesons cannot match with that of its quarks. 
However, this helicity suppression can be alleviated in the scalar meson production because of the 
non- vanishing orbital angular momentum L z with the scalar state. Consequently, weak annihilation 
contributions to B — > SP can be much larger than the B — > PP case. 

Finally, it is worth mentioning that we shall implicitly use the narrow width approximation in 
the calculation of the B decays into resonances; that is, we will neglect the finite width effect even 
for very broad resonances such as a and k states. Under the narrow width approximation, the 
resonant decay rate respects a simple factorization relation (see e.g. [49] ]) 

F(B -» SP^ PxP 2 P) = T{B -» SP)B(S -» P 1 P 2 ). (4.14) 

It has been shown in 49] that in practice, this factorization relation works reasonably well even for 
charmed meson decays as long as the two-body decay D — * SP is kinematically allowed and the 
resonance is narrow. The off resonance peak effect of the intermediate resonant state will become 
important only when D — > SP is kinematically barely or even not allowed. The factorization 
relation presumably works much better in B decays due to its large energy release. 



1. B -> / (980)K and a (980)K decays 

The decay mode B — * fo(980)K has been studied in [5fJ within the framework of the pQCD 
approach based on the Izt factorization theorem. It is found that the branching ratio is of order 
5 x 10~ 6 (see Fig. 2 of the second reference in [5(3]), which is smaller than the measured value by 
a factor of 3. 

The penguin-dominated B — > foK decay receives two different types of penguin contributions 
as depicted in Fig. ^ In the expression of B — ► foK decay amplitudes given in Eq. (HU), the 
superscript u of the form factor F Q reminds us that it is the u quark component of /o involved 
in the form factor transition [Fig. ^a)]. In contrast, the superscript s of the decay constant fj Q 
indicates that it is the strange quark content of /o responsible for the penguin contribution of Fig. 
I2[b). Note that and penguin terms contribute constructively to ir°K~ but destructively 
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(a) (b) 



FIG. 1: Penguin contributions to B~ -> f {980)K~. 




0° 90° 180° 

6 



FIG. 2: The branching ratio of B~ — > /o(980)-ftT _ versus the mixing angle of strange and non- 
strange components of /o(980), where the middle bold solid curve inside the allowed region corre- 
sponds to the central value. For simplicity, theoretical errors due to weak annihilation and hard 
spectator interactions are not taken into account. The horizontal band within the dashed lines 
shows the experimentally allowed region with one sigma error. 

to foK~ . Therefore, the contribution to B — > JqK from Fig. 1(a) will be severely suppressed. 
Likewise, the contribution from Fig. 1(b) is suppressed by f£° ~ mf /mi,- Hence, it is naively 
expected that the foK rate is smaller than the tt°K one. However, as shown in Appendix B, the 
scale dependent decay constant fi is much larger than f n owing to its scale dependence and the 
large radiative corrections to the quark loops in the OPE series. As a consequence, the branching 
ratio of B — ► foK turns out to be comparable to and even larger than B — » tt K. 

Based on the QCD factorization approach, we obtain B{B~ — * JqK") = (9.0 — 13.5) x 1CP 6 for 
25° < 9 < 40° and (12.0 - 17.2) x 10~ 6 for 140° < 6 < 165° (Fig. 2), where only the central values 
are quoted. 6 Hence, the short-distance contributions suffice to explain the observed large rates of 



6 The calculated branching ratios in the present work are slightly larger than that in |4£| because of the 
larger scalar decay constant /f and different estimates of the leading- twist LCDA for /o(980). It was 
originally argued in |4o| that while the extrinsic gluon contribution to B — > foK is negligible, the intrinsic 
gluon within the B meson may play an eminent role for the enhancement of fo(980)K. 
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(b) 

FIG. 3: Penguin contributions to B~ — ► /o(980)/T _ in the 4-quark picture for /o(980). 
f K- and f K°. 

Thus far we have discussed f$K modes with the two-quark assignment for the /o(980). It is 
natural to ask what will happen if fo is a four-quark bound state. Naively, one may wonder if the 
energetic /o(980) produced in B decays is dominated by the four-quark configuration as it requires 
to pick up two energetic quark-antiquark pairs to form a fast-moving light four-quark scalar meson. 
The Fock states of /o(980) consist of qq, q 2 q 2 , qqg, • • • , etc. It is thus expected that the distribution 
amplitude of fo would be smaller in the four-quark model than in the two-quark picture. Naively, 
the observed B — > fo(98Q)K rates seem to imply that the two-quark component of /o(980) play an 
essential role for this weak decay. 

Nevertheless, as pointed out in [j^], the number of the quark diagrams for the penguin contri- 
butions to B — > /o(980)/C (Fig. |3J) in the four-quark scheme for /q(980) is two times as many as 
that in the usual 2-quark picture (Fig. [J). That is, besides the factorizable diagrams in Fig. Ela), 
there exist two more nonfactorizable contributions depicted in Fig. |3Jb). Therefore, a priori there 
is no reason that the B — > fo(980)K rate will be suppressed if fo is a four-quark state. However, 
in practice, it is difficult to give quantitative predictions based on this scenario as the nonfactor- 
izable diagrams are usually not amenable. Moreover, even for the factorizable contributions, the 
calculation of the /o(980) decay constant and its form factors is beyond the conventional quark 
model, though an attempt has been made in |32j]. In order to make quantitative calculations for 
B — * fo(980)K, we have assumed the conventional 2-quark description of the light scalar mesons. 
However, as explained before, the fact that its rate can be accommodated in the 2-quark picture 
for /o(980) does not mean that the measurement of B — ► foK can be used to distinguish between 
the 2-quark and 4-quark assignment for /o(980). 

We next turn to B — ► ao(980)K decays. A main difference between o!qK and foK modes is 
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that the latter receives the dominant contribution from the s quark component of the /o [see Fig. 
mb)], while such a contribution vanishes in the former mode even when a[j is assigned with the 
ss(uu — dd)/V2 quark content. Because of the destructive interference between the 04 and a$ terms, 
the penguin contributions related to the u quark component of the ao and /o are largely suppressed. 
Consequently, the weak annihilation contribution becomes as important as the penguin one. For 
example, the branching ratio of a$K~ is of order 3.3 x 10~ 7 in the absence of weak annihilation, 
while it becomes 2.4 x 10 -6 when weak annihilation is turned on. From Table [V] we see that 
T(B —* clqK) <C T(B — > foK) and the OqK rate is enhanced by a factor of 2 for charged ao- The 
predicted central value of B(B° — ► a^K^) is larger than the current upper limit by a factor of 
2. However, one cannot conclude definitely at this stage that the 2-quark picture for ao(980) is 
ruled out since it is still consistent with experiment when theoretical uncertainties are taken into 
account. Nevertheless, as we shall see below, when the unknown parameter pa for weak annihilation 
is fixed to be of order 0.7 in order to accommodate the -K"g(1430)7r data, this in turn implies too 
large ao(980)K rates compared to experiment. There will be more about this when we discuss 
B ^(1430)vr decays. Note that the prediction of B(B~ % (980)iT°) = 15 x 1(T 6 made in 
[5ll | in the absence of the gluonic component is ruled out by experiment. 



2. B -> a (980)vr, / (980)vr decays 

The tree dominated decays B — > ao(980)-7r, /o(980)-7r are governed by the B — ► ao and B — > 
/q transition form factors, respectively. The fair rate is rather small because of the small uu 
component in the /o(980) and the destructive interference between 04 and 06 penguin terms. Since 
the B — * ao(980) form factor is predicted to be similar to that for B —* tt one according to the 
covariant light-front model (see Sec. III.C), it is interesting to compare B — * clqii decays with 
B — > 7T7T. First, B° —* ir + is highly suppressed. This means that the B° — B° interference plays 
no role in the cl^tt^ channels. Thus the decays B° — > ciQir^ are expected to be self-tagging; that 
is, the charge of the pion identifies the B flavor. Second, we see from Table [V] that the branc hing 
ratio B(B° -> a^O ~ 7.6 x 1(T 6 is slightly larger than B(B° -» vr+Tr") = (4.5 ± 0.4) x 10" 6 Q 
and that B(B~ -» vr-vr ) > B(B~ -» ajjvr-) > B{B~ -> a 7r°). 

Just as the a^K~ mode, the predicted branching ratio B(B° -> 0^(980)^^) = (8.2^;^;^^) x 
10 -6 exceeds the current experimental limit of 3.3 x 10 -6 by more than a factor of 2 (cf. Table IV|) . 
If the measured rate of Oqt:~ is at the level of (1 ~ 2) x 10~ 6 or even smaller, this will imply a 
substantially smaller B — > 00 form factor than the B — > tt one. Hence, the four-quark explanation 
of the ao (see Fig. @J) is preferred to account for the B — > ao form factor suppression. We shall 
see later that since ao(1450) can be described by the qq quark model, the study of ag~(1450)7r~ 
relative to a,t(980)7r~ can provide a more strong test on the quark content of ao(980). It has been 
claimed in [52| that the positive identification of B°/B° — ► a t (980)7r =F is an evidence against the 
four-quark assignment of ao(980) or else for breakdown of perturbative QCD. We disagree and we 
argue below that if the branching ratio of B — > ajj~(1450)7r - is measured at the level of 3 x 10~ 6 
and the aQ~(980)-7r~ rate is found to be smaller, say, of order (1 ~ 2) x 10~ 6 , it will be likely to 
imply a 2-quark nature for ao(1450) and a four-quark assignment for ao(980). 
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FIG. 4: Tree contribution to i? — ► a^(980)7r in the 4-quark picture for ao(980). 

In short, although it is unlikely that the penguin-dominated decay B — > /o-ftT can be used to 
distinguish between the 2-quark and 4-quark assignment for /q(980), the decays B — » aovr and ao-f^ 
may serve for the same purpose for ao(980). For example, the former mode is tree dominated and 
its amplitude is proportional to the form factor Fq" 10 which is suppressed in the four-quark model 



for do (980). It has been claimed in |17( that a best candidate to distinguish the nature of the ao 
scalar is B(B~ — > 7r°) as the prediction for a four-quark model is one order of magnitude smaller 
than for the two quark assignment. We see from Table IV1 that the branching ratio of this mode is 
only of order 2 x 10~ 7 even when ao(980) is treated as a 2-quark state. Experimentally, it would 
be extremely difficult to test the ao(980) nature from the study of Oq (980)7r°. 

It is commonly assumed that only the valence quarks of the initial and final state hadrons 
participate in the decays. Nevertheless, a real hadron in QCD language should be described by a 
set of Fock states for which each state has the same quantum number as the hadron. For example, 

|a+(980)> = rjud) + r u %\udg) + ^ZJudss) + .... (4.15) 

The possibility that ao(980) can be viewed as a bound state of four quarks at low energies, while 
its 2-quark component manifests at high energies is also allowed by current experiments. 

Note that the production of ao(980) in hadronic B decays has not been seen so far and only 
some limits have been set. In contrast, the oo(980) production in charm decays has been measured 
in several places, e.g. D° — » K°a{](980) and K~Oq (980) in the three-body decays D° — > K + K~K° 
[5^ . It is conceivable that the scalar resonance ao(980) in B decays will be seen at B factories 
soon. 



3. B -> £T *(1430)7r decays 

For weak decays involving scalar mesons above 1 GeV such as i£g(1430), ao(1450) and /o(1500) 
we consider two different scenarios to evaluate their decay constants and LCDAs based on the QCD 
sum rule method (see Appendices B and C): i) 1£q(1430), ao(1450), /o(1500) are treated as the first 
excited states of k, ao(980) and /o(980), respectively, and (ii) they are the lowest lying resonances 
and the corresponding first excited states lie in between (2.0 ~ 2.3) GeV. Scenario 2 corresponds to 
the case that light scalar mesons are four-quark bound states, while all scalar mesons are made of 
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two quarks in scenario 1. The resultant decay constants and LCDAs for the scalar mesons above 1 
GeV in these two different scenarios are summarized in Appendices B and C. The B — > 1Tq(1430) 
form factors in scenarios 1 and 2 can be found in Table llVl It should be stressed that the decay 
constants of Kq (1430), ao(1450), /o(1500) have the signs flipped from scenario 1 to scenario 2 as 
explained in footnote 9 in Appendix B. 

As mentioned in the Introduction, there exists a two-fold experimental ambiguity in extracting 
the branching ratio of B~~ — > i"Cg(1430) 7r~: Belle found two different solutions for its branching 
ratios from the fit to B + —* K + ir + ir~ events [3]. The larger solution is consistent with BaBar 



while the other one is smaller by a factor of 5 [see Eq. (|2,4j) ]. It appears that the larger of 
the two solutions, namely, B(B~ — > i^o(1430) 7r~) ~ 45 x 10~ 6 , is preferable as it is consistent 
with the BaBar measurement and supported by a phenomenological estimate in 18(. However, 
since B~ —* K°ir~ has a branching ratio of order 24 x 10 -6 [ijj, one may wonder why the K*qTT~ 
production is much more favorable than K°tt~ , while the K^ir mode is comparable to K°n (see 
Table HJ). 

— *o — o 

To proceed we consider the pure penguin decays B~ — ► K ir~ and B~~ — > K ir~ for the purpose 
of illustration. The dominant penguin amplitudes read [see also Eq. (|A6|) ] 

A(B- - i?*V) oc K " rpag)*^ fK*Ft(m 2 K*)(m 2 B ~ ml), 

A(B~ -> K°tt-) oc (<4 + r*al\ K f K F^ {m 2 K ){m 2 B - ml), (4.16) 

where we have neglected annihilation contributions for the time being. Although the decay constant 
of Xq(1430), which is 37 ± 4 MeV in scenario 2 [cf. Eq. (|B16|) ]. is much smaller than that of the 



kaon, it is compensated by the large ratio r* x ° = 8.9 at [/, = 2.1 GeV compared to r*f = 1.1. 



Since the penguin coefficient is the same for both Kq-tt and Kir modes, it is thus expected that 

T(Kq it~)/T{K 7t~) ~ 3.2 in the absence of the 04 contribution. When 04 is turned on, we notice 

— *o — 

that its contribution is destructive to K ir~ and constructive to K ir~ . In order to see the effect 

of 04 explicitly we give the numerical results for the relevant (ttKq) at the scale [i = 2.1 GeV 



and for o^(ttK) 



ai = 


1.417 -*0. 181, 


02 = 


0.673 -iO.lll, 


al = 


-0.199 - i0.009, 


a\ = 


-0.162 - i0.059, 


< = 


-0.0558 - zO.0163, 


a% = 


-0.0602 - i0.0039, 


a u 8 = 


(79.4 — i4.8) x 10~ 5 , 


a c 8 = 


(78.5 - i2A) x 10~ 5 , 


a u w = 


= (70 - »64) x 10~ 4 , 


a c w -- 


= (70 - i62) x 10~ 4 , 


ai = 


0.993 + i0.0288, 


a 2 = 


0.144- iO.lll, 


a\ = 


-0.0267 - i0.0183, 


a\ = 


-0.0343 - i0.0064, 


< = 


-0.0568 - i0.0163, 


«6 = 


-0.0612 - i0.0039, 




(74.4 -i4.5) x 10~ 5 , 


a c 8 = 


(73.6 — t'2.3) x 10~ 5 , 


a? = 


= (-208 + ^90) x 10~ 5 , 


a c w = 


= (-209 + m) x 10~ 5 



(4.17) 



(4.18) 

where scenario 2 has been used to evaluate both ^(ttKq) and (^(irK). Comparing Eq. (|4. 17f) and 
with Eq. (|4.18|) . it is evident that vertex and spectator interaction corrections to 01,02,04 and 010 
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for ttKq are quite large compared to the corresponding a^irK) due mainly to the different nature 

of the Kq LCDA. Note that the clq and as penguin terms remain intact as they do not receive 

vertex and hard spectator interaction contributions. Since the magnitude of cO^kKq) is increased 

— *o 

significantly, it is clear that the K tt~ rate eventually becomes slightly smaller tl 
to the large destructive contribution from o^(ttKq). Hence, we conclude that B(B~ 



1 x 10" 



B{B- 



K TT 



Kin) 



in the absence of weak annihilation contributions. 



^From Tables IV1 and IVTI it is clear that when weak annihilation is turned on, the K^tt rates are 
highly suppressed in scenario 1 due to the large destructive contributions from the defaulted weak 
annihilation. In order to accommodate the data, one has to take into account the power correc- 
tions due to the non-vanishing pa and pn from weak annihilation and hard spectator interactions, 
resepctively. Since power corrections are dominated by weak annihilation, a fit to the data yields 
PA ~ 0.4 for scenario 2 and pa ~ 0.7 for scenario 1, where we have taken <pA ~ 0. 



We see from Eq. ()A6j) that the amplitudes of B 



17*° - 



and B 



■o 



Kq 7r + are identical 



when the small contributions from the electroweak penguin and A u ai, X u b2 terms are neglected. 
This amounts to assuming the dominance of the AI = penguin contributions. Hence, these 
two modes should have the same rates under the isospin approximation |5J[. Likewise, T(B° — ► 
K* \ )/T(W -» K*--k+) = 1/2 is expected to hold in the isospin limit. Indeed, it is found in QCD 
factorization calculations that 7 



B(B° 



Ko°(1430)7T ) 



B(B° 
B(B~ 



-^0 



(1430)tt+ 
K^(U30)tt { 



B(B- 
t(B°)B{B- 



K* °(U30)tt- 



t(B-)B{B -> K*-(1430)tt- 



n ci +0.01+0.08+0.04 
u - oi -0.02-0.02-0.06 
+0.01+0.02+0.04 
-0.02-0.01-0.10 

+0.03+0.43+0.79 
-'-0.02-0.02-0.02 
n co+0.05+0.01+0.17 
u - oo -0.03-0.01-0.05 

n o 1 +0.06+0.62+0.93 
u - ol -0.07-0.04-0.57 
-0.05+0.03+0.18 
-0.06-0.03-0.27 



0.47 
0.30 



0.90 



scenario 1 
scenario 2 

scenario 1 
scenario 2 

scenario 1 
scenario 2. 



(4.19) 
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Consequently, the ambiguity in regard to B~ — * K ir~ found by Belle can be resolved by the 

Kq~tt + . As noted in passing, both BaBar and Belle measurements of 
B- -» K* u n- and B° -> Kq~ "7r + [see Eq. ()2.4|) and Table ILT| do respect the isospin relation. It 



measurement of B 



is also important to measure the ratio of B{W — > K q° (I430)ir°) / B(B U — > Kq~ (1430)7r + ) to see if 
it is close to one half. At any rate, both BaBar and Belle should measure all Kq(1430)tt modes 
with a careful Dalitz plot analysis of nonresonant contributions to three-body decays to avoid any 
possible ambiguities. 

We now turn to the implications of sizable weak annihilation characterized by the parameter pa 
which is of order 0.7 in scenario 1 and O(0.4) in scenario 2. We find that all the calculated ao(980)K 
rates are too large compared to experiment. For example, B(B° — ► a^(980)i^ _ ) ~ 31.4 x 10~ 6 for 
PA = 0.7 and ~ 14.6 x 10~ 6 for pa = 0.4. Both are ruled out by the current limit of 1.9 x 10 -6 . 



/-ftO 



7 From Table IVl and Eq. (|4.19l) . it appears that the mode K^~tt° does not respect the approximated isospin 
relation T(B~ — > K^tt )/T(B° — ► Kq~7t + ) = 1/2. This is mainly ascribed to the large cancellation 
between penguin and annihilation terms in the amplitude of B~ — > Kq~it [see Eq. (IA6(l ] and the 
remaining term proportional to (d2<5^ + 3(ag — 0,7) /2) breaks isospin symmetry. 
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This clearly indicates that oo(980) cannot be a purely two-quark state and that scenario 2 in which 
the light scalar meson is assigned to be a four-quark state is preferable. 



4. B -> a (U50)K, a (1450)vr decays 

For B — ► ao(1450)7r and an(1450)-fT decays, the calculated results should be reliable as the 
ao(1450) can be described by the qq quark model. Just as an(980)-K" modes, weak annihilation 
gives a dominant contribution to an(1450)if rates. It is found that their rates are much larger in 
scenario 1 than in scenario 2 due to the relative sign difference between the Gegenbauer moments 
B\ and B3 for ao(1450) and the sign of the ao(1450) decay constant flipped in these two scenarios 
(see Tables DO and [XT)) . The interference pattern between the penguin and annihilation amplitudes 
is generally opposite in scenarios 1 and 2. For example, the interference in 2?° — ► (1450)K~ is 
constructive in scenario 1 but becomes destructive in scenario 2. By the same token, the aQ~(1450)7r~ 
and OQ(1450)-7r° rates are also quite different in scenarios 1 and 2. 

As discussed in the previous subsection on ETg (1430)7r, predictions under scenario 2 are more 
preferable. Hence, if the branching ratio of B° — > a^(1450)-7r =F is measured at the level of 4 x 1CP 6 
and the ag~(980)7r~ rate is found to be smaller, say, of order (1 ~ 2) x 10~ 6 or even smaller than this, 
it will be likely to imply a 2-quark nature for an (1450) and a four-quark assignment for ao(980). 
Note that the naive estimate of 20 x 10~ 6 made by for this mode appears to be too large due 
to the usage of a large B — > ao(1450) form factor, ^^ a °( 1450 )(o) = 0.46. Experimentally, an (1450) 
will be more difficult to identify than an (980) because of its broad width, 265 ± 13 MeV J22j. 

5. B° — > k + K~ as spectroscope for k four quark state 

As for k (or ETq(800)), there is a nice and unique place where one can discriminate between 
the 4-quark and 2-quark pictures for the k meson, namely, the B° —* k + K~ decay. Recall that 
B° — > K + K~ is strongly suppressed as it can only proceed through the VF-exchange diagram. The 
experimental upper bound on its branching ratio is 0.6 x 10 -6 19J, |22] while it is predicted to be 



of order 1 x 10~ 8 theoretically (see e.g. j^]). Naively S° — > k + K~ is also rather suppressed if k is 
made of two quarks. However, if k has primarily a four-quark content, this decay can receive a tree 
contribution as depicted in Fig. E^b). Hence, if B° — > k + K~ is observed at the level of > 10~ 7 , 
it may imply a four-quark content for the k. Presumably, this can be checked from the Dalitz plot 
analysis of the three-body decay B° — > K + K~it° or B° — > K°K~ir + . As noticed before, scenario 
2 is more favored for explaining the B — > ^(1430)^ data. This already implies that k is preferred 
to be a four-quark state. 

Unlike the other light scalar mesons, the experimental evidence for k is still controversial. The k 
state has been reported by E791 in the analysis of D + — > K~tt + tt + with the mass 797±19±43 MeV 
and width 410 ±43 ±87 MeV 155J. However, CLEO did not see evidence for the k in D° -> R-tt+tt 
[5^. The k state was also reported by the reanalyses of LASS data on irK scattering phase shifts 
using the T-matrix method and the unitarization method combined with chiral symmetry [3 ■ 
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d u 

(a) (b) 

FIG. 5: Annihilation and tree contributions to 1?° — > k + K~ in the 4-quark picture for k. 

Most recently, BES has reported the evidence for the k in J/ip — ► X i^ + 7r~ process with the mass 
878 ± 23±H MeV and width 499 ± 52±|f MeV 0. 

It is interesting to notice that the decays B — > DfK~ and B° —* Z?* (2317) + i ; r~, the analogues 
of B° — > K + K~ and S° — ► n + K~, have been measured recently. The measured branching ratios 
are -D+if - ) = (3.8 ± 1.3) x 10~ 5 22] and £(5° -► D*+K-)B(D*+ -► D+vr ) = (5.3±^;g ± 



1.6) x 10 6 |60j. Since -D* (2317) + is dominated by the hadronic decay into Dfir , it is clear that 
T(W -> D*jK-) > T(W -► D+K-). These two decays can only proceed via a short-distance 
VF-exchange process or through the long-distance final-state rescattering processes 
DfK~ and 5° — ► Dq + 7t~ — > D*qK~. (In fact, the rescattering process has the same topology 
as ^-exchange.) Since B(lf -> D+tt') « 2.8 x 10~ 3 > B(3?° -» D+K~), it is thus expected 
that the decay i?° — > DfK~ is dominated by the long-distance rescattering process. As B(B° — ► 
D^+vr-) < 1.8x 10- 4 0, we will naively conclude that T(B° -> D*£K-)/T(I? DfK~) < 0.06, 
in contradiction to the experimental observation. Nevertheless, if -D* (2317) + is a bound state of 
csdd |olT |. then a tree diagram similar to Fig. Efb) will contribute and this may allow us to explain 
why T(B° -» L>: +if-) > r(5° - £>+#")■ 



6 1 . B ^ air decays 

The tree dominated -B — > (T7r decays are expected to have similar rates as B — ► 7r°7r ones if the 
ct meson is assumed to be a bound state of 2 quarks. Assuming that a has similar decay constant 
and LCDA as / (980), it is found that B(B~ cttt") ps 4.5 x 10" 6 and B(B° -> cttt ) ps 1.7 x 10~ 7 . 
The former is to be compared with the upper limit 4.1 x 10 -6 Q|. 



7. Direct CP asymmetries 

We see from Tables IVIll and IVlIII that CP partial rate asymmetries in those charmless B — > SP 
decays with branching ratios > 10~ 6 are in general at most a few percents. This is ascribed 
to the fact that the strong phases calculable in QCD factorization are generally small and that 
the observation of direct CP violation requires at least two different contributing amplitudes with 
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TABLE VII: Same as Table IVl except for CP asymmetries (in %). 



Mode 




Theory 


Expt 


Mode 




Theory 


Expt 


B- -> 


/ (980)iT- 


n ,.+0.0+0.0+0.6 
"•^-0.0-0.0-0.6 


-2 0+ 6 - 8 
z -"-6.5 


T?^ 


/o(980)^ U 


n 7+0.0+0.0+0.1 
"• ' -0.0-0.0-0.1 


-6 ±21 


B~ -» 


/ (980)vr- 


9 n+0.2+0.1+43.0 
z -"-0.2-3.4-41.9 


-50 ± 54 


B°^ 


/o(980)tt° 


qo C+H.8+ 2.9+ 44.1 
°°-"-ll. 0-18. 3-114.4 




B~ -> 


ag(980)if- 


q 0+1.6+1.5+52.8 
°-°-l.l-1.5-63.6 




Tf^ 


o^(980)A- 


q 4 +1.4+1.4+51.3 
°- -1.0-1.4-61.9 




B- -> 


a (980)1?° 


n Q+0.1+0. 1+1.2 
"• y -0.1-0. 2-0.9 




B°^ 


ao 1 (980)7?° 


n 7+0.1+0.1+1.0 
"•' -0.1-0.1-0.5 




5- -> 


ag(980)vr- 


n fl+0.1+0. 1+3.6 
"•"-0.1-0.2-3.8 




B° — > 


o^"(980)7r _ 


n q+0. 2+0.5+23.6 
"•°-0. 2-0.3-22.8 




J3~ -» 


a (980)vr° 


P.K n+6.1+7.8+24.6 
DO - y -8.3-5. 2-22.0 




U°^ 
B°^ 


a (980)vr+ 
a[j(980)vr° 


7K ^+ 5.9+4.0+21.0 
'"•°-10.5-5. 9-36.1 
oa q+12.3+ 9.1+28.6 
-J" 4 - 5 - 8.3-11.6-30.5 




B~ -» 


ag(1450)K- 


n Q+0.5+0.9+21.4 
"• y -0.3-0.6-18.7 




Tf^ 


0^(1450)^ 


n q+0.5+0.9+19.7 
"• J -0.3-0.6-18.9 




S~ -> 


a (1450)1?° 


n q+0.1+0. 1+0.2 

"••J-0.1-0. 1-6.6 






(1450)1?° 


n q+0. 1+0.1+0. 2 
"•°-0. 0-0.1-1. 7 




5- -> 


ag(1450)vr- 


9 q+0.2+0.4+5.4 
z - y -0.2-0.3-5.7 




B ^ 


0^(1450)^" 


n A +0.1+0.5+36.2 
"•^-0.1-0.3-35.4 




s- -> 


ao (1450)tt° 


1Q H+1.8+3.2+44.9 
iy -°-6. 4-3.6-46.6 




B° — >■ 
B°^ 


ao(1450)7r + 
ag(1450)vr° 


r Q 9+13.8+12.6+ 33.0 
oy - z -36. 7-66. 6-152. 7 
q 9 q+11. 7+7.8+66. 2 
oz '°-16. 7-6.6-52.4 




B- -> 


Zo°(1430)vr- 


a A+2.8+ 4.1+63.8 
^-4.9-25.7-30.5 


-5+ 5 

°-8 


B°^ 


Ko"(1430)7r+ 


1 r 1+5.0+18.5+16.1 
iJ - L -3.6- 5.9-21.2 


-7 ±14 


5- -► 


if^(1430)7r° 


A0 1+12.6+78.6+128.7 
^•"L— 10.1— 2.1- 12.3 




B°^ 


^ (1430)vr 


q a +0.4+0.3+10.8 
°-^-0.3-7.4- 9.1 


-34 ±19 



TABLE VIII: Same as Table fVTl except for CP asymmetries (in %). 



Mode 



Theory Expt 



Mode 



Theory Expt 



w 

B° 
B° 
B° 
B° 
B° 
B° 



B' 
B' 
B- 
B~ 

B~ 



ag(1450)K- 
a (1450)1?° 
ag(1450)vr- 
ao(1450)7r° 



i^o°(1430)^- 



ca 7+ 2.3+ 3.6+ 10.4 
-33.1-16.2-101.0 



B- -» i<: *-(1430)7r 



19.9 



4.1+ 8.8+8.4 
17.0-13.5-4.5 
n Q+0.2+0.1+6.3 
"• y -0.2-0.3-6.6 
a-, ,,+6.7+1.0+123.7 
^ LA -G.2-1.7- 64.8 



1.1 



+1.3+0.2+ 8.0 



0.7-0.1-17.8 
a 9+3.6+1.2+5.0 
^• z -2.9-1.2-3.7 



-,+5 



■ 0^(1450)^" 

■ a° (1450)1?° 

■ af(1450)7i- 

■ a (1450)vr+ 
• ag(1450)7r° 

■ Ko-(1430)vr+ 

■ l?o°(1430)vr 



10.0 

96.4 

0.6+0.5+1.9 
6.5-6.2-6.2 
q 7+0.0+0.3+46.9 



C-q q+ 4.8+ 5.3- 
°°-°-35.0-18.3 



4.9~ 



-0.0-0.2 
37 g+10.0+1.3 



45.2 

48.6 

11.9-2.1+68.0 
9H s+7.5+3.3+43.3 
z "-°-4.9-4.0-45.6 
q O+1.9+0.3+ 8.3 
°-°-3. 6-0.3-13.2 

n R +o. 5+0.1+0. 7 

"■"-0.3-0.1-0.7 



-7± 14 
-34 ± 19 



distinct strong and weak phases. Hence, if the observed direct CP asymmetry is of order 0(0.1) 
or larg er, then strong phases induced from power corrections could be important. As pointed out 
in 47| |. final-state rescattering processes can have important effects on the decay rates and their 
direct CP violation, especially for color-suppressed and penguin-dominated modes. However, this 
is beyond the scope of the present work. 



8. Mixing-induced CP asymmetries 

It is of great interest to measure the mixing-induced indirect CP asymmetries Sf for penguin- 
dominated modes and compare them to the one inferred from the charmonium mode (J /ipKg) in 
B decays. It is expected in the Standard Model that sin2/3 e fj defined via Sf = — rjf sin 2/3 e g with 
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TABLE IX: Mixing-induced CP parameter AS" = sin2/3 e g — sui2/3ckm m scenarios 1 and 2 as 
explained in Appendices B and C. The sources of theoretical errors are same as in previous tables 
except the last one is from the uncertainty in the unitarity angle 7. 



Mode 



Theory (Scenario 1) 



Theory (Scenario 2) 



B° -> / °(980)Ks 
W -> a%(980)K s 
B° -> a$(1450)K 3 



B° -> i^o°( 14 30)vr 



0.023 
0.022 
0.023 
0.004 



0.000+0.000+0.001+0.001 
0.000-0.000-0.001-0.001 

+0.000+0.000+0.005+0.001 
0.000-0.000-0.006-0.001 
+0.000+0.000+0.027+0.001 
0.000-0.000-0.001-0.001 
+0.005+0.010+0.030+0.000 
0.008-0.040-0.036-0.000 



n n01 +0.014+0.008+0.031+0.001 
U.UZJ._ .000-0.000-0.009-0.001 

n n9i+ - 001 + - 000 +°- 004 + - 001 

u - u ^ 1 -0. 002-0.013-0. 008-0.001 



rjf being the CP eigenvalue of the final state / should be equal to Sj/^k s with a small deviation at 
most O(0.1) - See j^] for recent studies of sin 2/3 e g in some of P — > PP and PU modes using the 
QCD factorization approach with or without the presence of final state interactions. In Table HXl 
we show the predictions on the mixing-induced CP parameter AS = sin 2/3 e ff — sin 2/3ckm for the 
CP eigenstates f (980)K s , a%(980)K s , a%(U5Q)K s and Ko°(1430)7r°, where only the CP compo- 
nent of i^Q°(1430) namely, Kstt , is considered in the last mode. In addition to the theoretical 
errors considered before, the uncertainty of 7° in the unitarity angle 7 is included. Note that main 
errors arise from the uncertainties in annihilation contributions and 7. Our results indicate that 
ASf in these penguin dominated modes are positive and very small. 

V. CONCLUSIONS 

In this work we have studied the hadronic B decays into a scalar meson and a pseudoscalar 
meson within the framework of QCD factorization. Vertex corrections, hard spectator interactions 
and weak annihilation contributions to the hadronic B — > SP decays are studied using the QCD 
factorization approach. Our main results are as follows: 

• Based on the QCD sum rule method, we have derived the leading-twist light-cone distribution 
amplitudes (LCDAs) of scalar mesons and their decay constants. It is found that the scalar 
decay constant is much larger than the previous estimates owing to its scale dependence and 
the large radiative corrections to the quark loops in the OPE series. Unlike the pseudoscalar 
or vector mesons, the scalar LCDAs are governed by the odd Gegenbauer polynomials. 

• While it is widely believed that light scalar mesons such as /o(980), ao(980), k are pre- 
dominately four-quark states, in practice it is difficult to make quantitative predictions on 
B — > SP based on the four-quark picture for S as it involves not only the form factors and 
decay constants that are beyond the conventional quark model but also additional nonfac- 
torizable contributions that are difficult to estimate. Hence, in practice we shall assume the 
two-quark scenario for light scalar mesons in calculations. 

• The short-distance approach suffices to explain the observed large rates of foK~ and foK° 
that receive major penguin contributions from the penguin process b — * sss and are governed 
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by the large fo scalar decay constant. When /o(980) is assigned as a four-quark bound state, 
there exist two times more diagrams contributing to B — > fo(980)K. Therefore, although 
the fo(980)K rates can be accommodated in the 2-quark picture for /o(980), it does not 
mean that the measurement of B — > JqK can be used to distinguish between the 2-quark 
and 4-quark assignment for /o(980). 

When ao(980) is treated as a qq bound state, it is found that the predicted 5° — > aQ~(980)7r~ 
and aQ~(980)i"C~ rates exceed substantially the current experimental limits. Hence, a four- 
quark assignment for ao(980) is favored. The ao(980)K and ao(1450)if receive dominant 
contributions from weak annihilation. 

Belle found two different solutions for the branching ratios of B + — ► Kq (1430)°7r + from the 
fit to B + — ► K + tt + tt~ events. The larger solution is consistent with BaBar while the other 
one is smaller by a factor of 5. Based on the isospin argument, we have shown that the 
smaller of the two solutions is ruled out by the measurements of Kq (1430)~7r + by BaBar and 
Belle. 

For B — > ao(1450)7r, ao(1450)K and A"q(1430)7t decays, we have explored two possible sce- 
narios for the scalar mesons above 1 GeV in the QCD sum rule method, depending on whether 
the light scalars re, ao(980) and /o(980) are treated as the lowest lying qq states or four-quark 
particles. We pointed out that in both scenarios, one needs sizable weak annihilation in order 
to accommodate the KqTt data. This in turn implies that all the predicted ao(980)K rates 
in scenario 1 will be too large compared to the current limits if do (980) is a bound state 
of two quarks. This means that the scenario in which the scalar mesons above 1 GeV are 
lowest lying qq scalar state and the light scalar mesons are four-quark states is preferable. 
The branching ratio of 5° — > a^(1450)7r T is predicted to be at the level of 4 x 10~ 6 . 

The decay a — ► k + K~ can be used to discriminate between the 4-quark and 2-quark nature 
for the re meson. This mode is strongly suppressed if re is made of two quarks as it can proceed 
through the VF-exchange process. However, if re is predominately a four-quark state, it will 
receive a color-allowed tree contribution. Hence, an observation of this channel at the level of 
> 1CP 7 would mostly imply a four-quark picture for the re. Presumably, this can be checked 
from the Dalitz plot analysis of three-body decay B — ► K + K~ir° or 

Direct CP asymmetries in those decay modes with branching ratios > 10~ 6 are usually 
small of order a few percents. However, final-state rescattering processes can have important 
impact on the decay rates and their direct CP violation. 

Mixing-induced CP asymmetries in the penguin dominated SP modes such as fo(980)Ks, 
a° Q (980)K s , 0^(1450)^5 and Zo (1430)[^:57r ]7r are studied. Their deviations from 
sin 2/?ckm are found to be positive (AS > 0) and tiny. 
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APPENDIX A: DECAY AMPLITUDES OF B —> SP 

The B — > SP (PS) decay amplitudes can be either evaluated directly or obtained readily 
from B — > VP (PV) amplitudes with the replacements: fv$v( x ) &s{%) an d my fy<& v (x) — ► 
— ms<& s c(x). (The factor of i will be taken care of by the factorizable amplitudes of B — ► SP shown 
below.) To make the replacements more transparent, it is convenient to employ the LCDA <&s( x ) 
in the form ()3.23j) and factor out the decay constants fs in $s{ x ) an d fs in [see Eq. ()3.25|) ]. 

so that we have 8 



*v(x) - $ s (z), $„(x) -» & s (x), fv^fs, »£--r|, (Al) 



where 



r v (n)- 2mv fvM , r s (a) - 2mS ^ - ^ (A2) 

x rn b (fi) f v ' x m h (n) fs m 6 (/x)(m 2 (/i) - mi(/i)) ' 

and use of Eq. Q3.1U[) has been made. For the neutral scalars a, fo and a(j, becomes divergent 
while fs vanishes. In this case one needs to express fs r x by fs^x w ^ n 

= (A3) 

With the above-mentioned replacements, the quantity Am 1 m 2 anci the coefficients of the flavor 
operators a P defined in 0] read 

Gp { {m 2 B - m 2 p )F BP (m 2 s )fs; for M X M 2 = PS, 

Mi M-z ^<y_ {m 2 B _ m 2 s)F BS {m 2 p)fp . {olMl M 2 = SP, 



o^{M x M 2 



a\{M x M 2 ) + a p 5 (MiM 2 ); for M X M 2 = PS, 
al (Mi M 2 ) - a p 5 {M 1 M 2 ); for Mi M 2 = SP, 



We found in the present work that it is most suitable to define the LCDAs of scalar mesons including 
decay constants. In this appendix we try to make connections between B — ► SP and B — > VP amplitudes. 
The latter have been worked out in detail in 42]. Since the LCDAs in |2| are defined with the decay 
constants being excluded, for our purposes it is more convenient to factor out the decay constants in the 
scalar LCDAs so that it is ready to obtain B — > SP amplitudes from B — > VP ones via the replacement 
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oP^M x M 2 ) 



c4 EW (M 1 M 2 ) 



f a? (Mi M 2 )-r^(M!M 2 ); for Mi M 2 = PS, 
\c%(M 1 M2)-r%c%(M 1 M 2 )\ for M l M 2 = SP, 

f ag(MiM 2 ) + a?(AfiM 2 ); for MiM 2 = PS, 
I (Mi M 2 ) - af (Mi M 2 ) ; for Mj M 2 = SP, 

f a? (MiM 2 ) - r^ag(MiM 2 ); for M X M 2 = PS, 
i < (MiM 2 ) -r^af(MiM 2 ); for M 1 M 2 = SP, 



(A4) 



where 



2m| 



' x m 6 ( / u)(m 2 + mi)(/i)' 

It should be stressed that the a p F^ p and a P F^ s terms in the decay amplitudes have an opposite 
sign. 

Applying the replacement (|A1|) and Eq. (|A4|) to the P> — ► yP and Py amplitudes given 
in Appendix of 42j, we obtain the following the factorizable amplitudes of the decays B — ► 
/ if, a 7r, <T7r, ao-^T, i^o^r 



p K p , v K p 
«4 - r X fl 6 + °10 - r X a S 



x f K F Bf ° (m 2 K )« ~ m%) + ( ag - Jag ) f/°j% ^ x (mJ,)(m|, - m&) 



1 



Kf s 



fB 



(h% + b 3 +b 3 



Gp 

71 



E A«{(4-r^-i(a? 

p=u,c ^ ^ 



10 'x 



x f K F Bf °(m 2 K )(m B -m%) + 

~ fB 



p p 

a£ - -a F o 



r f x ° 1% F BK (m%)(ml - ml) 



(h - ^h,Ew)fd K + (b 3 - 2 & 3,ew)k/ s 



A(B~ -» a° K' 



p K p , p K p 

a 4 — r x a 6 + a 10 — r x a 8 



A{B~ -» a iT 



= -TS^f(* 

p=-u,c ^ 

x feFo^-tmDtml - rojj - / B (fcdJ + 1,3 + f>3,Ew)„ *-J, 

V ^ p=u,c (. 



n p - r K n p - - (n v -r K n p ) 



x /kPo 500 ^)^! - m* ) - / B (fe 2 ^ + 63 + 6 3 ,Ew) aoi ,|, 
Gp 



^ E + < 

p=n,c ( 



r x a 6 + a 10 - r x a 



aoK 



x f K F Bao (m 2 K )(m 2 B - m 2 ao ) - f B (b 3 - -b 3 , E w) aoK 
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a°K°) 



G F 



E A H f«4-^6-^Ko-^)) 

f K F^(m 2 K ){m 2 B - m 2 J - f B (b 3 - h 3 , EW ) aoK 

V ^ p= n ,c I 

UF» fS (ml)(ml - m%) + fog - log) j% ^(mj,)^ - m*) 



/ovr ) 



13 13 

(6l<5g - 63 + -& 3 ,EW + 2 6 4,EW) + {blK - 63 + g 6 3,EW + -&4,Ew) u/ d 



- /b (M£ + h + h,Ew) f u n + (b 2 SZ + b 3 + h,Ew) nf u 

= % E 4 S) { ("^ + <- « - |K - «?) - - r;a§)) 

z p=u,c I v z z ' /(T 

x /^KX™! - m^) + (ag - ^ag) ff/jj, F^(m%)(m 2 B - m 2 w 
r 1 3 

+ /b 

--^W^-^-.j.e). 

V 1 p=u,c L ^O 71 " 

x (m*)(m| - r<) + (ag - ±ag) F ^«)(m| - m*) 

- /b (&2<S£ + fo 3 + ^ew)^ + {b 2 Sl + b 3 + 63,EW) W(T « 
= ^ E A J> S) { ("^ + < - rjag - |(ag - a?) - - rjag)) 

p=u,c L ^ ' 



cr vr ) 



x Aif ff6 (mJ)(m| - r<) + ( ag - Jag ) j f^ F ^(m^ )(m| - 



F°"0 f d P-Btt/ 2 



+ /B 



13 13 

(hS? - 63 + 2 6 3,ew + -64,Ew) a ^ + (h6 p a - 63 + t^ew + -^.ew)^ 



^ E ^ (a^ + ag-rjag + ^o-^ag) 

p= n ,c I 

x a ° (ml){m 2 B -m 2 J- f B [(b 3 +b A - i& 3 ,EW - ^4,Ew) a , )T 



ao7r 



a ?r 



+ (61^ + 64 + &4,EW) TO0 



^ E ^ + «4 " ^X°«6 + «?0 " tfag) 



7rao 



x f aa F^{ml){m 2 B - ml) + / B [(63 + 64 - ^3,EW - ^4,Ew) mi) 
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+ (&i<5£ + &4 + hjsw] 



A(B~ 



P IT V , V TV P 

a A — r x a 6 -+■ a 10 — r x a 8 



A(B- 



a TT 



0\ 



p=u,c { 

x UF^(ml){m 2 B -m 2 ao ) 

~ SB [(&2<5£ + ^3 + &3,Ew) ao7r " {b 2 SP + b 3 + 6 3 ,EW) 



ao7r 
2 



z / 7ra 



7rao 



p=«,c 

1 3 
+ /b f(&2<5£ + h + b3,Ew) wao - (W£ + 63 + &3,ew' 



ao 7r 



A(B° 



a°vr ) 



A(B- 



K Q 7T ) 



A(B~ - K *-7r ) 



= E W (<**S - aS + r^l + i« - rjag) + ^(a 9 - a 7 )) 

x /^(m^ml - m^) - (ag - lag) f«°/ ao (m| - ml)F^(m 2 J 

+ fB + 63 + 26 4 - ^3,EW " &4, EW )) ao . 

+ {hSP + 63 + 26 4 - ^(& 3 ,EW " &4,Ew)) ffao ] |, 

v 2 p=UjC { V 2 y 

x fK S Fj*(m 2 K .)(m 2 B - ml) + f B (b 2 5i + 63 + ^EW^j, 

= ^ E A W (*i«S + aS " r?ag + a? - r?a|) 

p=u,c ^ 

x /^F ^(m^)(m|-m2)- 
+ /b(62^ + 63 + 63,ew)^ 5 }. 

V ^ p=M, C I 

x fK*F^{m 2 K ,){m 2 B - ml) + / B (6 3 - ^3,ew)^*|, 



r x X + a io - r x °: 



A 7T ) 



p=u,c 



■kK* 
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x f K *Ft(rn 2 K »)(m' B -mi) 



+ /fl( - ^3 + 2 6 3,EW)^* >, 



A'* 



UF* K Hrnl)(m%-m 2 K *) 
(A6) 



where Ap 9 ' 1 = V pb V* q with q = d,s and 



2m^ 



m b (^)(m u (fi) + m s (/x)) ' 



rJ°(Ai) 



2m 



A ( r 



mb(ji)(m s (jj,) - rn q (fj,)) ' 



2m 



m b (n)(m d (p,) - m u {fj)) ' 



'(A*) 



2m, 



OQ 



m 6 (/i) 



2 m 



/o 



m b {(j)' 



(A7) 



Note that the /o~c mixing angle (i.e. sin#) and Clebsch-Gordon coefficient l/y/2 have been included 
in the /o(980) form factors F B ^o ' and decay constants fj^ d and likewise for the form factors F Ba o' 
and decay constants f% ,d - Throughout, the order of the arguments of the a^{M\M2) and bi{M\M2) 
coefficients is dictated by the subscript M1M2, where M2 is the emitted meson and M\ shares the 
same spectator quark with the B meson. For the annihilation diagram, M\ is referred to the one 
containing an antiquark from the weak vertex, while Mi contains a quark from the weak vertex. 



APPENDIX B: DETERMINATION OF THE SCALAR COUPLINGS OF 
SCALAR MESONS 

To determine the scalar decay constant fs of the scalar meson 5 defined by (0\q~2qi\S) = msfs, 
we consider the following two-point correlation function 

U(q 2 ) = i [ d 4 xe iqx {<d\T{f^{x)j q2qi \Q)\{)) , (Bl) 



with j q2qi = qiqi- The above correlation function can be calculated from the hadron and quark- 
gluon dynamical points of view, respectively. Therefore, the correlation function arising from the 
lowest-lying meson S can be approximately written as 

4ii = If^, ( B2) 

m S ~ 1 n •'O S — q z 

where n is the QCD operator-product-expansion (OPE) result at the quark-gluon level, so is 
the threshold of the higher resonant states, and the contributions originating from higher resonances 
are approximated by 

- / ds (B3) 

n J s s - q z 

We apply the Borel transformation to both sides of Eq. (|B2|) to improve the convergence of the 
OPE series and suppress the contributions from higher resonances. Consequently, the sum rule for 
lowest lying resonance with OPE series up to dimension 6 and 0(a s ) corrections reads jf^ 



m 2 p -rn%,M* ( ^ = 3 4 

SS \as(M)J 8tt 2 
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l,a s G 2 , fl \._ , fl ,. 

+ — > + I 2 mi +rri2 )( qiqi > + ( + mi i( q2q2 > 



^[\rn^q igs a ■ G qi ) + \ 



m 2 {qig s a ■ Gqi) + -mi(q 2 g s a ■ Gq 2 ) - ira s {q 1 a^ l ,\ a q 2 q 2 a flu X a qi) 



1 1 \ 

-Kas-iqu^q! q^X a qi) - Tra s -(q 2 ~f ll \ a q 2 q 2 ^\ a q 2 ) I , (B4) 

where /(l) = 1 - e~ s °/ M2 (1 + s /M 2 ), 1(1) = f^ so/M 2 (lni) ln(- ]nt)dt, the scale dependence of f s 
is 

f^=fM^mT' (B5) 



and the anomalous dimensions of relevant operators can be found in Ref. 65] to be 



a s Oo) \ b 



4 



{qq)» = (qq) 



Q! S (mo) \ 6 

" :1 1 n..;/,i J ' 



(5 s ga • = (^ s gcr • Gq)^ 



2_ 

:',b 



Qs(mo) V 

(a s G 2 ) M = (a s G 2 ) M , (B6) 

with b = (IIN C — 2nj)/3, where we have neglected the anomalous dimensions of the 4-quark 
operators. In the numerical analysis, we shall use a s (l GeV) = 0.517 corresponding to the world 
average a s (mz) = 0.1213, and the following values for vacuum condensates and quark masses at 
the scale /i = 1 GeV [65J: 

(a s G a ^G a ^) = 0.474 GeV 4 /(4^) , 

(uu) 9* (dd) = -(0.24 GeV) 3 , (la) = 0.8(uu) , 

(m u + m d ) /2 = 5 MeV , m s = 119 MeV , 

(g s uaGu) = (g s daGd) = — Q.8(uu), (g s saGs) = 0.8(g s uaGu) . 

We adopt the vacuum saturation approximation for describing the four-quark condensates, i.e., 

(0\qr i X a qqr i X a q\0) = --^Tr^F^T^ X a )(qq) 2 . (B8) 

Taking the logarithm of both sides of Eq. (|B4|) and then applying the differential operator M 4 d/dM 2 
to them, one can obtain the mass sum rule for the lowest-lying resonance S, where so is determined 
by the maximum stability of the sum rule. Substituting the obtained so and mass into Eq. (|B4|) . 
one arrives at the sum rule for the decay constant fs- 

Nevertheless, in order to extract the decay constant fs> for the first excited state S' , we shall 
consider two lowest lying states on the left hand side of Eq. (jB4|) . i.e., 

m 2 sp s e~ m s/ M2 + m 2 ,p s , <r™\>l M * = - f° ds e"^ Imn OPE ( S ) . (B9) 

vr Jo 
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1. a (980) and a (1450) 



Taking qiq2 = ud and considering only the ground state meson, we obtain 
m s k (0.99 ± 0.05) GeV, 

f s {l GeV) = 370 MeV, f s (2.1 GeV) = 440 MeV , (B10) 

corresponding to sq ~ 3.1 GeV 2 and the Borel window 1.1 GeV 2 < M 2 < 1.6 GeV 2 , so that 
the resulting mass is consistent with ao(980). However, if one would like to have the mass re- 
sult of the ground state to be consistent with that of ao(1450), then one should choose a larger 
so — 6.0 GeV 2 together with the Borel window with a larger magnitude: 2.6 GeV 2 < M 2 < 
3.1 GeV 2 . Since k, ao(980) and /o(980) may be four-quark states, we therefore explore two 
possible scenarios: (i) In scenario 1, we treat re, oo(980), /o(980) as the lowest lying states, and 
Kq (1430), ao(1450), /o(1500) as the corresponding first excited states, respectively, where we have 
assumed that /o(980) and _/o(1500) are dominated by the ss component and (ii) we assume in 
scenario 2 that Kq (1430), ao (1450), fo (1500) are the lowest lying resonances and the corresponding 
first excited states lie between (2.0 ~ 2.3) GeV. Scenario 2 corresponds to the case that light scalar 
mesons are four-quark bound states, while all scalar mesons are made of two quarks in scenario 1. 

In the numerical analysis, we adopt the first two lowest resonances as inputs in these two 
scenarios and perform the quadratic fits to both the left-hand side and right-hand side of the 
renormalization-improved sum rules in Eq. ()B9|) . We find that in scenario 1 the resulting threshold 
and Borel window are so = (5.0±0.3) GeV 2 and 1.1 GeV 2 < M 2 < 1.6 GeV 2 , respectively, while in 
scenario 2, s = (9.0 ± 1.0) GeV 2 and 2.6 GeV 2 < M 2 < 3.1 GeV 2 . Thus for a (980) and a (1450), 
we obtain 

/«o(980)(l GeV) = (365 ± 20) MeV, / oo(B80) (2.1 GeV) = (450 ± 25) MeV, 
/«o(l460)(l GeV) = -(280 ± 30) MeV, f M14BO) (2.1 GeV) = -(345 ± 35) MeV, 

in scenario 1 and 

/«o(l460) (1 GeV) = (460 ± 50) MeV, / ao(1450) (2.1 GeV) = (570 ± 60) MeV, 
fs>(l GeV) = (390 ± 80) MeV, f s ,(2.l GeV) = (480 ± 100) MeV, 

in scenario 2, where S' denotes the first excited state. Note that the sign of the decay constants 
for the excited states in scenario 1 cannot be determined in the QCD sum rule approach [see Eqs. 
(B.9) and (C.6)]. They are fixed from the signs of the form factors as shown in Table llVl using the 
potential model calculation. 9 



9 In the quark model with a simple harmonic like potential, the wave functions for a state with the quantum 
numbers (n, I, m) is given by fniip 2 / P 2 )Yim{P) exp(— p 2 /2/3 2 ) up to an overall sign, with fio(x) = fu(x) = 
1 and f2i[x) = y/b/2(l — 2x/5). For the n = 2, 1 = 1 state, the decay constant fs is dominated by the 
second term in /21, while the B — > S form factors is governed by the first term in /21 as the spectator 
light quark in the B meson is soft. Consequently, the decay constant and the form factor for the excited 
state have opposite signs. The overall sign with the wave function can be fixed by the sign of the form 
factor which is chosen to be positive in general practice. 
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2. / (980) and / (1500) 



Here we will assume that /o(980) and /o(1500) are both dominated by the ss component, i.e. 
j ss = ss. The results read 

/ /o(980) (l GeV) = (370 ± 20) MeV, // ( 98 o)(2.1 GeV) = (460 ± 25) MeV, 

4(1500) (1 GeV) = -(255 ± 30) MeV, / /o(1500 ) (2.1 GeV) = -(315 ± 35) MeV, 
in scenario 1, and 

4(1500) (1 GeV) = (490 ± 50)MeV, 7/ (l500)(2.1 GeV) = (605 ± 60) MeV, 

f s >(l GeV) = (375 ± 80) MeV, f s >(2.1 GeV) = (465 ± 100) MeV, 

in scenario 2. 



3. k(800) and £T *(1430) 

The relevant current is j qs = qs with q = u or d for the cases of k(800) and -?Tq(1430). Using the 
single resonance approximation as given in Eq. (|B4|) . we find that the lowest lying mass roughly 
equals to (0.86 ± 0.02) GeV 2 , corresponding to so — 2.4 GeV and the Borel window of 0.8 GeV 2 < 
M 2 < 1.3 GeV 2 . In analogy with the case of ao(1450), if i^o(1430) is justified by the result of 
the lowest lying mass sum rule, then it is necessary to have a large threshold so — 6.0 GeV 2 
corresponding to a larger Borel mass region 2.6 GeV 2 < M 2 < 3.1 GeV 2 , where the stable plateau 
can be reached. 

For k(800) and ^(1430), we find 

/ K(800) (l GeV) = (340 ± 20) MeV, 7«( 8 oo)(2.1 GeV) = (420 ± 25) MeV, 
/tf 5 (i430) (1 GeV) = -(300 ± 30) MeV, 7tf*(l430)O = 2-1 GeV) = -(370 ± 35) MeV, 
in scenario 1 and 

7^(1430) (1 GeV) = (445 ± 50) MeV, 7**(1430) (2-1 GeV) = (550 ± 60) MeV, 

f S '(l GeV) = -(420 ± 80) MeV, /s'( 2 -l GeV) = -(520 ± 100) MeV, ' ' 

in scenario 2. 

Two remarks are in order. First, if neglecting the RG improvement for the mass sum rules 
and considering only the lowest lying resonance state, the results, as stressed in Ref. j64|, become 
sensitive to the values of the four quark condensates for which the vacuum saturation approximation 
has been applied. Then it is possible to have results to be consistent ao ( 1450), 1430) and 
/o(1500) in the range of 0.8 GeV 2 < M 2 < 1.2 GeV 2 if so is larger and four-quark condensates are 
several times larger than that in the vacuum saturation approximation. Second, thus far we have 
considered renormalization-group (RG) improved QCD sum rules. It is found that sum rule results 
become insensitive to four-quark condensates if the RG improved effects are considered. For the 
RG improved mass sum rules, if taking ao(1450), Kq (1430) and /o(1500) as lowest resonances, then 
it is necessary to have a large threshold sq ^ 4.9 GeV 2 corresponding to a much larger Borel mass 



region 2.6 GeV 2 < M 2 < 3.2 GeV 2 , in contrast with the conclusion in Ref. 36] where the stable 
Borel window for the JCq(1430) mass sum rule is 1.0 GeV 2 < M 2 < 1.2 GeV 2 . 
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APPENDIX C: LEADING TWIST LCDAS FOR SCALAR MESONS 



The LCDA $s(x, /i) corresponding to the quark content q\q 2 is defined by 

Jo 



(01) 



where x (x = 1 — x) is the momentum fraction carried by the quark q (antiquark q) and \i is the 
normalization scale of the LCDA. $5(3:,//) can be expanded in a series of Gegenbauer polynomi- 
als 0,0 



$sO,m) = 7s6x(l - x) 



(C2) 



^^(M)Cf /2 (2x-l) 
L/=o 

where multiplicatively renormalizable coefficients (or the so-called Gegenbauer moments) are given 
by 



(C3) 



which vanish for even I in the SU(3) limit. Consider the following two-point correlation function 



d 4 xe^(0|T(O;(x)Ot( )|0) = (zq) l+l k(q 2 ) , 



(04) 



(2x - l) l <f> s (x)dx = (zp) l+1 f s {e s ) 



(05) 



H(?) = « 

where 

<0|0,|^(P)> = <0|«2 ftiz D) l qi\S(p)) = (zp) l+1 
(Q\O\S(p)) = (0\q 2 q 1 \S(p))=m s f s , 

with z 2 = and £ = 2x — 1. 

We shall saturate the physical spectrum with two lowest lying resonances for reasons to be 
explained later. Therefore, the correlation function I\ can be approximately written as 

' ' ' 2 vrio 



■;??, 



5' 



So WpPE/ a ) 

as — ^ — 



(06) 



where S 1 and S' refer to the lowest and first excited resonance states, respectively, and 

f\ (<Z2<?2> , 10/ - 3 (q 2 g s a ■ Gq 2 



hii 2 ) 



m 00 + m 



91 1 m q2 m qi 



+ 



16tt 2 ^ 1 + 2 I + 1 

1(41-5) (g 2 G 2 )(q 2 q 2 



In 



18 



+ (-1) 



z+i 



16tH 



+ 



m q2 + m, 
1 + 2 



qi 



24 

m, 



</2 



m, 



91 



Z + l 



In 



(gigi) 10/ - 3 {qig s a • Ggi) 

o "T 



1(41-5) (g 2 s G 2 )(q iqi ) 



q* 24 q A 18 

In terms of the above defined moments the sum rule reads 



(07) 



3 M 2 | 



16tt 2 



m q2 + m. 
1 + 2 



qi 



+ 



m 
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m, 



l + l 



/(0) + (g 2 g 2 ) + 



10/ - 3 (q 2 g s a ■ Gq 2 ) 



24 



| 1(41-5) (g 2 s G 2 )(q 2 q 2 ) + , +1 



36 



+<9i9i) + 



M 4 

10/ - 3 (qig s a ■ Gqi) 



</2 



M 2 



24 



M 2 



+ 



16tt 2 "" ^ Z + 2 
Z(4/-5) (g 2 G 2 )(q iqi )' 



l + l 



1/(0) 
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with /(O) = 1 — e s °/ m2 ; while the Gegenbauer moments are given by 



*"<"> - ^w^ lc ' 2{( ° )y (C9) 



Conformal invariance in QCD indicates that partial waves in the expansion of Ht?s(x, /J,) in Eq. ()C2}) 
with different conformal spin cannot mix under renormalization to the leading-order accuracy. 
Consequently, the Gegenbauer moments B\ renormalize multiplicatively: 



B l (ri=B l ( f i )[-^) , (CIO) 



where the one-loop anomalous dimensions are |6 

^-^('- (FFW+q +^^l- (C11) 

with Cf = (N 2 — 1)/(2N C ). Note that JsBq is independent of the renormalization scale. It should be 
also stressed that if only the lowest resonances are taken into account in Eq. QC8|) , the resultant mass 
reading from the sum rule that follows the same line as before by taking [(M 4 <9 / 'dM 2 ) In] to both 
sides of Eq. (|C8|) is less than 0.4 GeV, which is too small compared with the observables. Therefore, 
in the numerical analysis, we shall consider the first two lowest resonances and perform the quadratic 
fits to both the left-hand side and right-hand side of the renormalization-improved moment sum 
rules, given in Eq. (|C8l) . within the Borel window M min < M 2 < M max with M^ in ,M max G 
(1.1 GeV 2 , 1.6 GeV 2 ) [and M min ,M max G (0.8 GeV 2 , 1.3 GeV 2 )] corresponding to (£ ao ,/ ) [and 
(£k,A"cJ(143o)] m scenario 1 and M min ,M max G (2.6 GeV 2 , 3.1 GeV 2 ) in scenario 2, where the Borel 
windows are same as those in the previous section. It should be noted that for the moment sum 
rule for (£ ) in the large I limit, the actual expansion parameter is M 2 /l. Therefore, for (£ 3 ) we 
rescale the Borel windows to be M min , M max G (1.4 GeV 2 , 1.9 GeV 2 ) for a , f [and M min , M max G 
(1.1 GeV 2 , 1.6 GeV 2 ) for k, ifg(1430)] in scenario 1 and M^ in ,M max G (2.9 GeV 2 , 3.4 GeV 2 ) in 
scenario 2. Furthermore, for I > 5 and fixed M 2 , the OPE series are convergent slowly or even 
divergent, i.e. the resulting sum-rule result becomes less reliable. Following the same line as given 
in the previous section, we explore two possible scenarios. The results for the fist and second 
moments of together with the fist and second Gegenbauer moments are collected in Tables IXl 
and EH 
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TABLE X: Gegenbauer moments at the scales fj, = 1 GeV and 2.1 GeV (shown in parentheses) in 
scenario 1. 



State 


\s/ \s / 


Bi B-i 


a (980) 


-0.56 ±0.05 -0.21 ±0.03 


-0.93 ±0.10(-0.64± 0.07) 0.14 ±0.08 (0.08 ± 0.04) 


a (1450) 


0.53 ±0.20 0.00 ±0.04 


0.89 ±0.20 (0.62 ±0.14) -1.38 ± 0.18(-0.81 ± 0.11) 


/o(980) 


-0.47 ±0.05 -0.20 ±0.03 


-0.78 ±0.08(-0.54 ± 0.06) 0.02 ± 0.07 (0.01 ± 0.04) 


/o(1500) 


0.48 ±0.24 -0.05 ±0.04 


0.80 ±0.40 (0.47 ±0.28) -1.32 ± 0.14(-0.77 ± 0.08) 


k(800) 


-0.55 ±0.07 -0.21 ±0.05 


-0.92 ±0.11(-0.64± 0.08) 0.15 ±0.09 (0.09 ± 0.05) 


i^ *(1430) 


0.35 ±0.07 -0.08 ±0.06 


0.58 ±0.07 (0.39 ±0.05) -1.20 ± 0.08(-0.70 ± 0.05) 



TABLE XI: Same as Table IXl except for scenario 2. 



State 




Bi B 3 


a (1450) 


-0.35 ±0.07 -0.24 ±0.06 


-0.58 ± 0.12(-0.40 ± 0.08) -0.49 ± 0.15(-0.29 ± 0.09) 


higher resonance 


0.44 ±0.27 0.22 ±0.11 


0.73 ±0.45 (0.51 ±0.26) 0.17 ±0.20 (0.10 ±0.12) 


/ (1500) 


-0.29 ±0.06 -0.19 ±0.05 


-0.48 ± 0.11(-0.33 ± 0.08) -0.37 ± 0.20(-0.22 ± 0.12) 


higher resonance 


0.34 ±0.30 0.16 ±0.15 


0.56 ±0.50 (0.39 ±0.35) 0.07 ± 0.23 (0.04 ±0.13) 


#5(1430) 


-0.35 ±0.08 -0.23 ±0.06 


-0.57 ± 0.13(-0.39 ± 0.09) -0.42 ± 0.22(-0.25 ± 0.13) 


higher resonance 


0.25 ±0.11 0.12 ±0.05 


0.41 ±0.34 (0.28 ±0.24) 0.09 ± 0.14 (0.05 ± 0.08) 
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